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Whereas preheating after chaotic and hybrid inflation models has been abundantly studied in the 
literature, preheating in small field inflation models, where the curvature of the infiaton potential 
is negative during inflation, remains less explored. In these models, a tachyonic instability at the 
end of inflation leads to a succession of exponentially large increases and decreases of the inflaton 
fluctuations as the inflaton condensate oscillates around the minimum of its potential. The net effect 
is a competition between low-momentum modes which grow and decrease significantly, and modes 
with higher momenta which grow less but also decrease less. We develop an analytical description of 
this process, which is analogous to the quantum mechanical problem of tunneling through a volcano- 
'n— «' shaped potential. Depending on the parameters, preheating may be so efficient that it completes in 

less than one oscillation of the inffaton condensate. Preheating after small field infiation may also be 
5— ( ' followed by a long matter-dominated stage before the universe thermalizes, depending on the energy 

scale of inflation and the details of the inflaton interactions. Finally, another feature of these models 
. is that the spectrum of the inflaton fluctuations at the end of preheating may be peaked around the 

Hubble scale. In fact, because preheating starts when the second slow-roll parameter \ri\ becomes of 
order unity while the flrst slow-roll parameter e is still much smaller than one, the universe is still 
inflating during preheating and the modes amplified by the initial tachyonic instability leave the 
Hubble radius. This may lead to an abundant production of primordial black holes and gravitational 
waves with frequencies today which are naturally small enough to fall into the range accessible by 
high-sensitivity interferometric experiments. 

QJ- PACS numbers: 98.80.Cq, 98.70.Vc 

^ ■ I. INTRODUCTION 

! Simple inflationary models can be broadly but conveniently classified into three main categories, see e.g. [l|, ^ 
•/"J I for reviews. One category comprises hybrid inflation models [3|, where the inflaton is responsible for the slow-roll 
dynamics while another scalar field is responsible for the end of inflation. Interesting aspects of these models include 
the possible production of topological defects at the end of inflation and possible realisations in supergravity as F- 
04 term or D-term inflation [1] and in string theory in the framework of brane inflation [2] . Another category comprising 
^""^ , chaotic inflation ^] includes the large field inflation models. They have the simplest field theoretic origin as a single 
■ massive scalar field is enough to generate infiation. Another relevant feature of these models is the large, super- 
. . I Planckian value of the infiaton and therefore the possibility to generate a significant amount of primordial gravity 
^ • waves during infiation. Unfortunately, such large values of the infiaton may render these models more difficult to 
realize in supergravity or string theory. In this respect, the category of small field or hilltop infiation models 0, 
which includes for instance the original new infiation model 6], may be more natural. In these models, inflation 
occurs near a maximum or an inflection point of the potential, where the curvature is negative. This makes the 
slow-roll conditions somewhat easier to achieve and leads to cosmological perturbations with a negative spectral index 
as favoured by observations. Another interesting aspect of these models is that inflation can occur at very low energy 
scales, while still generating an acceptable spectrum of primordial density perturbations 0- These models offer also 
string theoretic realisations, for instance in the form of racetrack models (4I where the imaginary part of a Kahler 
modulus acts as the inflaton. Another consequence of small field infiation follows from the Lyth bound on the infiaton 
excursion which implies that gravitational waves generated during infiation are highly suppressed in this setting. 

At the end of infiation, the infiaton condensate must decay and reheat the universe. It has become clear during the 
last twenty years or so ^ that, in most models, reheating starts with an explosive and non-perturbative production of 
large, non-thermal fiuctuations of the infiaton and other bosonic fields coupled to it, in the process of preheating [To| . 
The subsequent dynamics are characterized by a highly non-linear and turbulent-like evolution, before the system 
eventually settles into thermal equilibrium. Preheating may have many interesting consequences in cosmology, like 
the production of stochastic backgrounds of gravitational waves flTI- [T5| , primordial black holes p^l - fTgj , non-gaussian 
curvature perturbations [20I [2l| or primordial magnetic fields |22l |. Preheating after chaotic and hybrid inflation 
has been abundantly studied in the literature, both analytically and with numerical lattice simulations [23l - [25j . In 
models of chaotic inflation, preheating occurs typically via broad parametric resonance [Toj, where fields coupled to 
the inflaton are produced non-adiabatically as the inflaton oscillates around the minimum of its potential, or via 
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tachyonic resonance [26*1 in the presence of trilinear or non-renormalizable bosonic interactions for the inflaton. In 
models of hybrid inflation, field fiuctuations are amplified by a tachyonic (spinodal) instability [13] as the fields roll 
towards the true minimum of the potential at the end of inflation. 

By contrast, preheating after small field inflation remains less explored. On the analytical side, the problem is 
complicated by the fact that closed-form solutions for the evolution of the homogeneous background during preheating 
are in general not available, while standard approximations like WKB for the evolution of the perturbations are not 
applicable in this context. On the numerical side, lattice simulations of preheating in these models must cover scales 
that range from the Hubble rate at the end of inflation up to the inflaton mass at the minimum of the potential, 
which typically differ by several orders of magnitude. Preheating after new inflation was studied in [1^^, where it 
was described as a combination of both tachyonic amplification and non-adiabatic resonance. In this sense, it is 
somehow intermediate between preheating after chaotic inflation and preheating after hybrid inflation. Some aspects 
of preheating in this model where also studied numerically in [19] and in the context of Kahler moduli / Roulette 
inflation in [29|. The combination of tachyonic and non-adiabatic effects occurs also in preheating via trilinear or 
non-renormalizable interactions after chaotic inflation [2^. However, we will see that preheating after small field 
inflation occurs in a qualitatively different way. 

The main purpose of this paper is to develop the analytical understanding of preheating in the class of small field 
inflation models. This is a first necessary step before a further study of the dynamics with lattice simulations and 
of the cosmological consequences of preheating in these models. To illustrate some of the qualitative features of 
preheating after small field infiation, consider the inflaton potential shown in Fig. [T] Slow-roll inflation ends when 
the second slow-roll parameter \r]\ becomes of order one, where the curvature of the potential is negative and of the 
order of the Hubble rate squared. The inflaton then rolls down the potential, passes through an inflection point 
where the curvature vanishes and oscillates around the minimum with an amplitude controlled by Hubble friction. 
Before the inflaton reaches the inflection point, its fluctuations have a negative effective mass squared and modes with 
momentum < —V"((j)) are amplified by a tachyonic instability. This process may affect a very wide range of scales, 
from the Hubble rate at the end of inflation up to the maximum of ]V"(^)], but it is much more efficient for the low- 
momentum modes. The tachyonic amplification stops when the inflaton crosses the inflexion point of the potential, 
oscillates around the minimum and goes back towards the tachyonic region. During this interval of time, the mass 
of the inflaton varies very non-adiabatically with time. This can lead to a further growth of the perturbations, but 
we will see that this effect is typically negligible. However, the evolution of the modes during this time interval when 
the inflaton oscillates around the minimum of the potential has dramatic consequences for the fate of the fluctuations 
when the inflaton goes back into the region where the curvature is negative. Indeed, during this second tachyonic 
episode, the amplitude of the fluctuations starts to decrease exponentially with time. For the modes with sufficiently 
low momenta, this exponential decrease occurs during the amount of time the infiaton condensate climbs back along 
its potential. If the inflaton climbed back up to the point from where it started, this effect would exactly compensate 
the growth of these modes during the first tachyonic episode. By contrast, the modes with higher momentum, which 
were amplified at a slower rate during the first tachyonic episode, continue to be amplified during the second one. 
The whole process may then repeat itself several times as the infiaton oscillates around the minimum of the potential. 
The net effect then follows from a competition between low-momentum modes which grow a lot but also decrease a 
lot and modes with higher momenta which grow less but also decrease less. As we will see, the analytical description 
of this process is very similar to the problem of tunneling through a volcano-shaped potential (two potential barriers 
separated by a deep potential well) in quantum mechanics. Alternatively, the initial tachyonic growth may be so 
efficient that preheating ends during the first tachyonic episode, depending on the parameters. 

As we will see, another specificity of small field models is that the spectrum of the infiaton fiuctuations amplified 
at the beginning of preheating is peaked around the Hubble rate, because the modes with k/a < H are amplified 
from the beginning of preheating when V" ^ —H^ while the modes with higher momenta are amplified only later and 
at a slower rate. The situation is different for preheating after chaotic inflation, where the fluctuations amplified by 
parametric resonance are typically peaked at scales k/(aH) > g^/"* 3> 1 with the initial resonance parameter q 3> 1 
for preheating to be efficient flO] . Similarly, for preheating after hybrid inflation, the typical scale of the fluctuations 
amplifled by the tachyonic instability is set by the curvature of the potential in the direction of the symmetry- 
breaking field and must be sub-Hubble to satisfy the so-called waterfall condition 0]. In fact, in models of small 
field infiation, the modes amplified by the initial tachyonic instability leave the Hubble radius, because preheating 
starts when the the second slow- roll parameter ]77] becomes of order unity while the first slow- roll parameter e is still 
much smaller than one, so that the universe is still inflating at the beginning of preheating. When preheating ends 
in less than one oscillation of the inflaton condensate, large density perturbations of the Hubble size may lead to an 
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FIG. 1: Schematic form of the inflaton potential in a model of small field inflation. 



abundant production of primordial black holes, see also |19l |. which may put strong constraints on models of small 
field inflation. This is in contrast with preheating after chaotic and hybrid inflation, where the field inhomogeneities 
are peaked at scales that are much shorter than the Hubble radius [IMllli. Furthermore, large field fluctuations at the 
Hubble scale and the fact that small field inflation can occur at very low energy scales may lead to the production of 
gravitational waves from preheating with frequencies today that are sufficiently small to fall into the range accessible 
by high-sensitivity ground-based and even space-based interferometric experiments 13|. A detailed study of these 
cosmological consequences of preheating after small field inflation will appear elsewhere [30| . 

The rest of the paper is organized as follows. In section HH we specify the models of small fleld inflation that we will 
consider and we describe the inflationary dynamics and the subsequent evolution of the inflaton condensate. Section 
mil is dedicated to the first stage of preheating, when the infiaton condensate first rolls towards the minimum of the 
potential. In Section ITVl we study the second stage of preheating, when the infiaton condensate oscillates around the 
minimum. Some of the results used in that Section are derived in the Appendix. In Section |Vl we first consider other 
possible non-perturbative decay channels for the condensate. We then discuss the perturbative decay of the infiaton 
after preheating and the resulting reheat temperature. Finally, Section I VII contains a summary of our results and 
directions for future work. 



II. SMALL FIELD INFLATION AND BACKGROUND EVOLUTION 



We consider small field inflation models with a single inflaton fleld (j). Inflation occurs near a maximum or an 
inflection point of the potential, which can be chosen as the origin (f> — 0. Denoting by M the energy scale during 
inflation, the potential can be parametrized as 



1-/ 



(1) 



where /(O) = /'(O) = 0. Restricting to positive values of cj), inflation occurs as the inflaton grows away from the 
origin. We denote hy (f) — v the minimum of the potential where 1^ = 0,/(1) = 1 and /'(I) = 0. The higher order 
derivatives of f{x) at x = 1 are usually of order one. In small fleld models, inflation occurs when (p ^ v ^. The 



If V > mpi, inflation may end when the inflaton is relatively close to the minimum of the potential <p v. This occurs for instance in 
modular inflation |31| and natural inflation [33l . These models belong also to the hilltop category but they are not small-field models 
(the infiaton takes instead Planckian values) and we do not consider them here. Preheating in these models is more similar to preheating 
after chaotic infiation [lC| . 
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inflationary potential may then be expanded into a power series around the origin as 



/(-)=- (-Y + - for <P^0 (2) 



P 

with p > 1 and A > 0. The absence of the p = 1 term in ^ may follow naturally from the fact that (j> — could 
be a fixed point of some internal symmetries, but this is not always the case. The location of any fixed point of 
internal symmetries is not important for the inflationary dynamics but it has crucial consequences for reheating, as 
we will see in Section |Vl The original new inflation model @ is well described by the potential ([1] [2]) with p — 4 
(up to logarithmic corrections) and </> = a fixed point of internal symmetries. In the original scenario, the inflaton 
was supposed to be a Higgs field in Grand Unified Theories with v ^ 10^^ GeV, but this is not necessary. Another 
example is provided by the model of MSSM inflation of [s^, where the inflaton corresponds to a flat direction of the 
Minimal Supersymmetric Standard Model. This model is well described by the potential ([l] [2|) with p — 3 and cj) — v 

a fixed point of internal symmetries ^, where v ^ (mA'^)^^^ is related to the cutoff scale A and to the soft SUSY 
breaking scale m ^ TeV. 

The details of preheating will of course depend on the explicit expression of the function f{4>/v) and in particular 
on the form of the potential around its minimum. However, in this section, we will first review the inflationary 
dynamics, in order to determine the initial conditions for preheating and the constraint to impose on M and v from 
the normalization of the CMB fiuctuations. For this, it is sufficient to consider the form ([1] [2|) of the potential in the 
vicinity of = 0. Let us first suppose that a single power of the infiaton in ^ dominates the potential throughout 
inflation. The slow-roll parameters then read 
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where a prime denotes derivative with respect to and mpi = I/VSttG is the reduced Planck mass. It follows that 
e/\ri\ ~ \?]\ (/)^/m|[, so that e <C |?7| during slow-roll (|?7| < 1) for small field models {(p mp\). The slow-roll phase 
then ends when \r]\ 1, while inflation continues afterwards until e ^ 1. Since jr/l must increase with (f> for the slow-roll 
phase to end, Eq. (|4]) implies that the inflaton potential must be dominated by a power (j>P with p > 2 at the end of 
slow-roll. This does not forbid the presence of terms with smaller powers of the infiaton in the potential (jll2p during 
inflation, but these terms must be subdominant at the end of slow-roll. In particular, a mass term 0^ is usually present 
in the potential when the model is embedded into supergravity and its coefficient must be somewhat smaller than its 
generic value in order to have jr/j < 1 during inflation, the so-called ry-problem It is possible that a mass term 0^ 
dominates the potential when cosmological scales leave the Hubble radius while a higher order term dominates at the 
end of inflation. However, it turns out that the constraint on M and v that follows from the normalization of CMB 
fluctuations in this case is practically the same as the one obtained when a single power c/)'^ with p > 2 dominates the 
potential throughout inflation, Eq. ([9|) below. Therefore, it is sufficient to consider a single power of the inflaton <j>P 
with p > 2 in the potential (jll2l) in order to determine both the initial conditions for preheating and the constraint 
to impose on M and v in order to obtain an acceptable spectrum of density perturbations. 

Preheating will start at the end of the slow-roll phase, i.e. when 77 ~ —1 but e <C 1, because the potential term 
then starts to dominate in the equation of motion of the infiaton fiuctuations, see the following Section. Eq. (|4]) gives 



1 

Hp - 1) m-pi 



(5) 



for the value 0e of the infiaton at the end of slow-roll when ~ —1. Using 3iJ(/) ~ — y leads to 
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/r2 / „, \P/(p-2) 
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(6) 



^ In I33II , the inflaton is normalized to </) = at the minimum of the potential. 
* The r;-problem in the context of small-field inflation is discussed in e.g. 0, Is^ . 
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for the time derivative of the infiaton at that time. These expressions provide the initial conditions for the inflaton 
condensate at the beginning of preheating. However, inflation will continue afterwards because e ^ 1 when 77 ~ —1. 
Inflation ends when e ~ 1, so that Eq. (jS]) gives 



A mpi 



(7) 



for the value 0^ of the inflaton at the end of inflation. Note that 0e ^ 4'e for v ^ mp\. When the inflaton condensate 
rolls from <j) = (pe to (f> = (f>e, preheating occurs while the universe is still inflating. 
The value (pN of the inflaton fleld N e- folds before the end of inflation is given by 

^)'"^A(p-2)A.4i (8) 

which can be derived from N — rrip^ J^^ d(j)/\/2e where the integral is dominated by its lower limit. The normalization 
of the CMB fluctuations requires 



V 



1/4 



~ 6.6 lO^*' GeV EE Minfl (9) 

when the cosmological scales leave the Hubble radius at iV = A^*. Using Eqs. (I3l8p . this gives 

/^yp-4)/(p-2) [(p_2)Ar^]-(^^i)/(^-^) Mg,a .... 
\mpj - y2Ai/(p-2) ■ ^ > 

For instance, for p = 3 and 4, this implies that M <C w. Depending on the values of v and p, inflation can occur in 
a very wide range of energy scales. Note also that Eqs. (jU and ((T0| hold in the usual scenario where the curvature 
perturbation is generated from the vacuum fluctuations of the inflaton field. If another field participates in the 
production the observed CMB fluctuations, like for instance in the curvaton scenario [s^, then the "~" should be 
replaced by a " <" in these two equations. Finally, deflning 

C EE V2 Ai/(f-2) _ 2)N,]^P-''^/<-P-^^ (11) 
the CMB normalization ([TU]) can be re-written 

^^c^r^^ypfJL) . (12) 

mpi \ Alinfl / 

For A ^ 1, Eq. ([TT]) gives C > 1. Eq. ([T2]) then implies that M < Mjnfl for any p> 2 and v < mpi. 

If the curvature perturbation is generated by the quantum fluctuations of the inflaton fleld only and if the term ^ 
with p > 2 dominates when cosmological scales leave the Hubble radius, the spectral index ti^ ~ 1 + 2?] — 6e is given 
by 

In small fleld models, inflation may occur at low energy scales and thermalization after preheating may be delayed 
(see Section |V|, so that the number of e- folds N^, before the end of inflation when cosmological scales leave the 
Hubble radius may be relatively small. This in turn may be constrained by the observational lower bound on 1 — 
via Eq. (|13p. However, contrary to Eq. (|10p . Eq. ([T^ is significantly modified if a mass term in 0^ dominates the 
potential when cosmological scales leave the Hubble radius, and the spectral index cannot be used to constrain the 
energy scale of inflation in that case. 

Examples of small field inflation models are described by the effective potential 



which reduces to (|1I2I) with X — 2p during inflation. In [7] , models of this kind with 

yP w M« m^7* (15) 
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were constructed in supergravity, where p > 2 and q are integers that depend on the symmetries of the theory. Using 
([T5|) in ([TU]) shows that, even in the usual scenario for the origin of the cosmological perturbations, inflation may occur 
at a very low energy scale M. The examples considered in [3| include a case where M corresponds to the intermediate 
scale of SUSY breaking in gravity-mediated scenarios, M ~ 10^^ GeV with v ^ 10^^ GeV, and a case with M ~ w ~ 1 
GeV. 

Let us now discuss the evolution of the inflaton condensate (t){t) starting from the end of slow-roll. It satisfies the 
Klein-Gordon equation 

(^-f 3F0 + !/'((/)) =0 (16) 

where a dot denotes the derivative with respect to the cosmological time i, a prime denotes the derivative with respect 
to (f> and H = d/a is the Hubble rate. A quantity of prime interest is the curvature of the potential V"{4>), which 
governs both the time scale for the variation of the inflaton condensate and the typical scales that will be amplified 
by preheating. At the beginning of preheating (/) = (/ig, 77 ~ — 1 so that 

l/"(0e)^-3i/f (17) 

where H,, is the Hubble rate at the end of slow-roll. At the minimum of the potential, (/) ~ v, we have instead 

V"{v) « ^ (18) 

where we have taken |/"(1)| sa 1. There must be an inflection point where V" = in between 0e and v, and we will 
denote it hy (j> = (j>i. Since jy" ((/>)| is increasing at the end of slow-roll, there is also a value of in between 0e and 
(pi where |y'(0)| is maximum, and we will denote it as ^ = Thus 

U"'(0„O = O and U"((/),) = (19) 

with < (pm < 4>i < V. For /"(I) ~ /'"(I) ~ 1, we have (pm 4>i ^ "v, whereas (j)e v for v <s:_ mpi. Note also 
that \V"{(j)e)\/V"{v) ~ v'^/nipi <C 1, so the inflaton stays a long time close to 0e before plunging rapidly towards the 
minimum of the potential. 

The condensate then oscillates around the minimum with an amplitude that decreases with time due to Hubble 
friction. The decrease can be estimated as in [2^ , from the conservation of energy that follows from Eq. ((T| 



i 02 + v{^) = ^ 0? + VicP,) - AV (20) 

where the dissipation term 

AU = 3 / dtH<p^ (21) 



results from the Hubble friction. Denoting by (pj the value of (f> after j complete oscillations of the inflaton (i.e. when 
is minimum and (pj = 0, while pj=Q — (pe), we have 



V{^,-i)~Vi^j)^3 / d^H^. (22) 

The integral above is dominated by its contribution around (p ^ v and can be approximated by taking H ~ 
Af2/(\/3mpi), (j) « and Acp sa 2v. During the first oscillations of the inflaton, the potential at (pj <C v can 
be approximated by (jll2p . Eq. (1^^ then gives 



V 



2^3 



A 




p V 



(23) 



after j > 1 inflaton oscillations. 
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III. FIRST STAGE OF PREHEATING: INFLATON'S ROLLING 



In this section, we study the amphfication of the inflaton fluctuations during the first stage of preheating, when the 
inflaton condensate first rolls towards the minimum of the potential at the end of slow-roll, i.e. when 4>e 1^ 4> 
We then estimate the conditions under which preheating ends during this first stage. The second stage of preheating 
(if preheating does not end during the first stage), when the inflaton condensate passes the minimum of the potential 
and oscillates around it, occurs in a qualitatively different way. This second stage of preheating, and the subsequent 
decay of the inflaton fluctuations into different forms of particles, will be considered in the following sections. 

The first stage of preheating depends essentially on the form of the potential at the end of inflation, so we can 
consider the generic form ((Tl [2]) for small field inflation without specifying the form of the potential around its 
minimum. The typical scale of the fluctuations amplified during this first stage will be of the order of the Hubble 
scale. In this case, the coupling of the inflaton fluctuations 5(f) to metric perturbations may be non-negligible. During 
the linear stage of preheating, we can consider the Mukhanov variable 

^ = J0 + $ (24) 

H 

where $ is the metric perturbation in the longitudinal gauge 

ds^ = - (f + 2$) dt^ + a^it) (1 - 2*) (25) 

with $ = ^ (since the anisotropic stress vanishes at linear order in the perturbations). In Fourier space, it satisfies 
the equation 

i^k + 3Hijk+ (^^ + V"i(b)+2^-2^ + 6H^ iPk^O (26) 

where k denotes the comoving wave-number. During inflation, we have the usual amplification of perturbations by 
the quasi-exponential expansion of the scale factor 



This provides the initial conditions for the evolution of the perturbation during preheating. 
Defining as usual v — c?l^ ip, Eq. (pS)) reduces to the equation of an harmonic oscillator 

Vk+ujl{t)vk^O (28) 

with frequency squared 

^fcW = ^ + ^"('/')+A (29) 

where A ~ H^. Preheating starts when the potential term V"{(j)) starts to dominate over A in this equation, i.e. at 
the end of the slow-roll phase when </) = t/fg, 77 ~ — 1 and V" ~ —37?^. Since V" < 0, the modes with k^/a^ < \V"\ 
have a negative frequency squared and they are amplified by a tachyonic instability. During all the time that the 
condensate rolls from = 0e to = ^ 0e, see Eq. ([7]), preheating occurs while the first slow- roll parameter e is 
still much smaller than one and the universe is still inflating. As a consequence, H is still approximately constant 
throughout the first stage of preheating. Eq. then reduces to the equation of motion for the infiaton fluctuations 
and the metric perturbation is negligible, ip ~ 54>. Furthermore, since the universe is still inflating, some of the modes 
that are amplified by the tachyonic instability leave the Hubble radius at the beginning of preheating. This is specific 
to small field infiation, where preheating starts when V" ~ —H^ and e <C f . 

As the inflaton condensate rolls towards the minimum of the potential, all the modes with k"^ /a? smaller than 
I V"((/))| are amplified by the tachyonic effect. This varies from |F"(0e)| ~ H"^ at the end of inflation, up to |l^"((/)m)| 
(the maximum of |F"(0)|, see dH])) which is typically of the order of |^"(0m)| V"{v) - Af^/w^ > ij2^ see dH]). For 
V <C ?7ipi, this represents a very wide range of momenta. However, the process is dominated by the low-momentum 
modes, which are amplified during a much longer period of time and at a much higher rate. Indeed, the maximal 
growth occurs for the modes with fc^/a^ < \V"\ ~ iH^ at the end of slow-roll, which are amplified from the beginning 
of preheating. In the following, we denote by — M^/(A/3mpi) the Hubble rate at the end of slow-roll / beginning 
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of preheating(t = te) and we normalize the scale factor as a{te) = 1. In this section, we will focus on the modes with 
k <^ He, which dominates the first stage of preheating. The evolution of the modes with higher momenta, which we 
will see dominate the second stage of preheating, will be studied in the next section. 

We can study the evolution of the low- momentum modes as in [2^ , by noticing that they behave as the velocity of 
the homogeneous inflaton condensate which is given by the conservation of energy (|20|) . Indeed, the time derivative 
of Eq. gives 



y + 3Hy+{V"i^)+3Hj y = (30) 

where y — (f). This reduces to Eq. ([26)) for fc^ <J ^ | V" \ (except at the very beginning of preheating and very close 
to the inflection point where V" vanishes), f/jfe + V" 'ipk — 0. The general solution for -0^, with k < then reads 

/t 
02 ^^^^ 

where ci and C2 are constants. As the inflaton rolls from cj) — (j)e to (j) — v, (/) and the integral above is a slowly 
varying function of time, so that we have essentially ipk oc </>. Note that this implies in particular that, for H constant, 
the comoving curvature perturbation TZ — Hip/cjjocHis conserved for super-Hubble modes, as it should. Indeed, for 
H constant and k/a ^ 0, Eqs. ([26l) and (|30|) are identical and ip oc (j> exactly. 

The initial velocity of inflaton condensate at the beginning of preheating, t — t^, is given in ([5]). Let us denote by 
t ~ ty the time when the inflaton condensate passes for the first time by the minimum of the potential (j) — v. Eq. (|20p 
gives <P{ty) ~ V2M^. Therefore 

t44^m^VE[Xip- ir-^] '"'-''^ (rj^Y"'"' (32) 

V'fe(ie) 0(te) L ^ J \ y ) V ; 

for the low- momentum modes. For v <^ mpi, the growth of the perturbation from t — te to t — ty is very large. Using 
Eq.dail) for Vfe(ie), Eq. ^ gives 

k^H.ity)\^^3[xip-ir^r^^^-^^ (r^)2./(.-.) j.^^^^ ^^^^ 

for k <, He- We have solved the equations of motion (fTHl [26]) for the inflaton condensate and the perturbations 
numerically for the potential (fM)) with p = 4, starting from the initial conditions at the end of slow-roll. The 
spectrum of k^ \ipk{ty)\'^ for different values of v is shown in Fig. [21 For the low-momentum modes k <^ H^, which 
dominate the spectrum, the results agree very well with the analytical prediction ()33[) . The amplitude of the spectrum 
is maximum at /c* ~ 4 He and then quickly decreases for larger values of k (this part of the spectrum will be studied 
in the next section). Because aH is still growing at the beginning of preheating, the peak of the spectrum at the end 
of the first stage is in fact slightly outside the Hubble radius, fc*/(a7?) ~ 0.3. 

Eq. p3|) is valid in the linear regime of preheating, when the amplitude of the inflaton fluctuations is much smaller 
than the inflaton condensate, |(50(t„)| ^ v. In that case, the energy density at t = ty is still dominated by the inflaton 
condensate, which starts oscillating around the minimum of the potential. Otherwise, the backreaction of the inflaton 
fluctuations shuts off the tachyonic effect and preheating ends in less than one oscillation of the condensate. We 
can estimate the variance of the inflaton fluctuations S4> ^ i/j s^t t = ty from Eq. ([55[) . (Scj)^) cx J dink k^ IV'feP- As 
usual, the integral must be regularized or restricted to a given range of scales due to the mild logarithmic divergence 
of the inflationary perturbations in the IR (in addition of course to the usual power law divergence of the vacuum 
fluctuations in the far UV). The integral is then dominated by the peak of the spectrum at k ^ He- This gives 

{Sc^'ity)) « [k'\Mtr)\'],^^^ « [Xip~ir-'f''-'' Hi (!!^)^^/(^-^) (34) 

where we have neglected factors of order unity. Preheating does not end during the flrst stage if {6(j)'^{ty)) < w^, i.e. 

[A (p-l)-r^^^^"''- • (35) 



In that case, the subsequent oscillations of the inflaton condensate lead to a second stage of preheating, as discussed 
in the next section. 
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FIG. 2: Spectra (^:^^ 



^ fc'' |?/>(f^)|^ as a function of k/He in the linear regime of preheating for the model (I14p with p = 4 



and t; = 10 GeV, 10 GeV and 10 GeV. The 3 spectra lie on top of each other and are in very good aggreement with 
Eq. (133). 

On the other hand, if the condition ([35)) is not satisfied, preheating ends already during the first stage. In that case, 
the spectrum of the inflaton fluctuations at the end of preheating is peaked around the Hubble scale, or even slightly 
outside the Hubble radius, before the energy gets redistributed by rescattering. As mentioned in the introduction, 
large field fluctuations at the Hubble scale are rather specific to preheating after small field infiation and may have 
interesting consequences for the production of primordial black holes and gravitational waves. Let us now discuss the 
kind of models where this occurs, i.e. for which the condition ([55]) is violated and preheating ends in less than one 
oscillation of the condensate. In the model of New Inflation [6^, p = 4, AI/v ~ 10~^ and A ^ 1, so preheating ends in 
less than one oscillation for v/mpi ^ 10~^ [2^. In the model of MSSM inflation of [s^l, p = 3, A ~ 1, M y/m v and 

V ^ (toA^)^^"' where m is a soft SUSY breaking mass and A the cutoff scale (denoted as A"^"^ = A/p '^/A„ in 33]). 
In that case, the condition psp is violated, and therefore the spectrum of inflaton fluctuations is peaked around the 
Hubble scale at the end of preheating, for A/mpi ^ (m/mpi)^^^. For m ^ TeV, this occurs for A/mp\ <; 10~^, in 
particular for A of the order of the GUT scale. 

In general, we can also determine the condition for preheating to end during the first stage in terms of v or M only 
by using their relation ([T0| from the CMB normalization. Eq. ([SS]) is then violated for 



Using (|12|) . this upper bound on v gives an upper bound on A/, so that preheating ends during the first stage if 
inflation occurs at a sufficiently low energy scale. This upper bound on M is very sensitive to the value of p, see 
Section|V]for details. In terms of v, the condition (p6)) for preheating to end in less than one oscillation of the inflaton 
condensate is roughly v/nipi < 10~^ — 10^^. 



When the condition ([55)) is satisfled, the inflaton condensate still carries most of the energy density at the end of 
the first stage of preheating and it starts to oscillate around the minimum of the potential. In this section, we study 
the evolution of the inflaton fluctuations during this second stage of preheating when v/mp\ > 10"^. The decay of 
the inflaton into other fields will be studied in the next section. 

In order to study the second stage of preheating quantitatively, we need to specify the form of the potential around 
its minimum. We will first consider the potential (1141) with p — A 




(36) 



IV. SECOND STAGE OF PREHEATING: INFLATON OSCILLATIONS 




(37) 



The generalization to other models of small field inflation will be discussed in sub-section IIVDI For the model p7p , 
the value 0e of the condensate at the end of slow- roll, ([S|), the infiection point and the point (j>m where —V" is 
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maximum, see (jl9[) . are given by 



1/4 



and 



1/4 



/24mpi ' \7J " V7. 

Note again that <^ (j)m ^ 4'i '^v. The mass of the inflaton at the minimum of the potential reads 

^ ]ur2 

and the maximum momentum of the inflaton fluctuations amplified by the tachyonic effect is 

4 



= V-V"i^n^) = 3777 ~ 0.435m; 



(38) 



(39) 



(40) 



To study the evolution of the inflaton fluctuations 
number 



71/4 V 

! "0, it is convenient to consider their comoving occupation 



Vk\ 



, .2 



(41) 



which is conserved when the frequency u)k{t) varies adiabatically with time (and is real). Here the modes Vk satisfy the 
oscillator equation (pS)) with the frequency (|29p . Although the occupation number is only well-defined when > 0, it 
is convenient to generalize the definition (PT|) to the case where < too, by using an "effective frequency squared" 
given by where V" is replaced by its absolute value. This reduces to the correct expression for the occupation 
number when w| > and it allows to also follow the evolution of the modes when they are tachyonic. 
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FIG. 3: Evolution with time of the occupation number 
njj for the mode with k = He (blue) and of the square 
of the inflaton condensate 4''^{t)/v'^ (red), for the model 
H37[) with expansion of the universe and v = 10^" GeV. 
Here Uk has been rescaled by an arbitrary factor in order 
to fit on the plot. Note that the plot does not display the 
evolution from the end of inflation but focus instead on 
the first oscillations of the condensate. 



FIG. 4: Evolution with time of the occupation number 
Uk for the modes with k = 10'^ He (blue), k = W'^ He 
(red) and k ^ 5x10^ He (black), for the model ((ST]) with 
expansion of the universe and v — 10^^ GeV. 



We solved numerically Eq. ([25]) for the model (1571) during the linear stage of preheating for different values of k and 
v/mpi, and we computed the occupation number as defined above. In Fig. [3l we show the evolution with time of nk 
for k = He and of the square of the inflaton condensate (jP{t)/v'^, during the first oscillations of the condensate, for 
V = 10^® GeV. The condensate first rolls slowly from (j) = (p^ v and then oscillates rapidly around the minimum 
(j) = V. During each oscillation, varies between a maximum value cj) ~ 2^/*w and a minimum value 0j which agrees 
well with and which increases at each oscillation due to Hubble friction. The main observation from Fig. [3] is 
that, for the modes with k ^ He which dominate the first stage of preheating, while the occupation number grows 
exponentially by a large amount when the condensate first rolls towards (f) = v, it also decreases exponentially by a 
large amount when the condensate rolls back towards the origin. In other words, the amplification of this mode during 
the first tachyonic episode is partly canceled by the second tachyonic episode. It would be fully canceled if during the 
second tachyonic episode the infiaton condensate rolled back up to the value from where it started. The process then 
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repeats itself as the condensate rolls back towards the minimum and oscillates. Overall, the occupation number of 
the modes with k/a <^ H oscillates with a large amplitude during the second stage of preheating, with almost no net 
growth after the first stage. Indeed, these low-momentum modes still follow the behaviour of the condensate velocity, 
Vk oc (j) and Wfc oc </> oc V'{(j)). When is maximum, the first term in (|41l) vanishes but the second one is maximum, so 
that Hk is maximum. When is minimum, the first term vanishes and the second one is small, so that nk is minimum. 

In Fig. [3l we show the evolution with time of the occupation number for modes with higher momenta. All these 
modes have k < fcmax in the tachyonic regime and are amplified by the tachyonic instability. As we discussed in the 
previous section, the modes with higher momenta grow less when the condensate first rolls towards the minimum of 
the potential. However, we see in Fig. [3] that these modes also decrease less when the condensate goes back towards 
the origin, so they progressively catch up with the modes with smaller momenta. Overall, there is a competition 
between low-momentum modes which grow a lot but also decrease a lot and modes with higher momenta which grow 
less but also decrease less ^. As the condensate oscillates around the minimum of the potential, the spectrum of 
the inflaton fluctuations is peaked at much larger momenta than during the first stage of preheating and the peak is 
progressively shifted towards the UV. Such a UV peak was attributed in [ll] to a non-adiabatic production of the 
inflaton fluctuations, but we will see that it follows rather from the large exponential decrease of the low-momentum 
modes when their frequency squared is negative. 

We will call tachyonic oscillation this succession of exponential increases and decreases of the low-momentum modes 
with k < /cmax- In order to study this process in more detail, let us first neglect the expansion of the universe, a{t) = 1. 
In that case, the inflaton condensate starts from some initial value (j)o <SC w, with an initial velocity that we can take 
to vanish, — 0. It then evolves periodically in time and comes back to its initial value (f>o after each oscillation. It 
is convenient to rescale the field, the time coordinate and the wave-numbers as 

, i^^t , ~k=^k. (42) 

V 2 mmin 

The conservation of energy ([20)1 with (po = Q then reads 

i 02 + v{4>) = V{M (43) 
where a dot now denotes the derivative with respect to the resettled time t and we have defined the rescaled potential 

V^l[l-^^y . (44) 



In terms of the rescaled variables (|42|) . the mode equation for the inflaton fluctuations (1281) without expansion of the 
universe reads 

ijk + Cblit) Vk ^ with Cblit) = P + V" (j){t)j (45) 
where a prime now denotes the derivative with respect to (fi. Defining 

~ko = \/-V"{M-^^o, (46) 

the modes with k < ko have a negative frequency squared initially and are amplified by the tachyonic instability from 
the beginning. On the other hand, the maximum wave-number amplified by the tachyonic instability ()40p is 



fcmax = (^y j ^ 0.87 . (47) 

In the rest of this section, we will always work with the rescaled variables HH HSl) : but from now on we omit the 
tilde on (j>, t, k and ujk to simplify notations. 

The mode equation (|45p can be viewed as a Schrodinger equation in one spatial dimension t, where the time evolution 
of the modes Vk{t) corresponds to the spatial profile of the wave- functions of states with energy k^ in the effective 



^ A similar process has been observed numerically for a potential V = Vo — 0^ -(- A in the context of preheating after hybrid inflation 
in the second paper of |27|| . 
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FIG. 5: The shape of the effective Schrodinger potential —V" [(j>{t)) in Eq. (|45|) for the model (1441) during one oscillation of 
the inflaton condensate. Here for convenience the coordinate t has been normalized to f = at the center of the "crater" - 
i.e. when the inflaton condensate reaches its maximum value before oscillating back. The full potential is made of a periodic 
succession of these volcano potentials during each oscillation of the inflaton condensate. 



potential ~V" {4>{t))- For the model (|44l) . the shape of this effective Schrodinger potential during one oscillation of 
the inflaton condensate is shown in Fig. [S] The full potential is made of a periodic succession of these potentials 
during each oscillation of the inflaton condensate. We see from Fig. [5] that the Schrodinger potential has a volcano 
shape, with two "crests" separated by a deep "crater". The time evolution of the modes Vk{t) during one oscillation 
of the inflaton condensate corresponds to a tunneling through this potential for k < fcmax and to a scattering above 
this potential for k > fcmax- In this language, the large exponential decrease of the low-momentum modes that we 
observed above when the inflaton condensate rolls back towards small values of 4> corresponds to a resonant tunneling 
of low-energy states ^. Indeed, the wave- function of these low-energy states is localized at the center of the deep 
"crater" and decreases symmetrically away from it. 

Let us first consider the time evolution of the modes Vk(t) during a single oscillation of the inflaton condensate, 
starting from <j) — <j)o ^ 1. We will normalize the time coordinate such that t = corresponds to the time after half a 
period of oscillation, when the inflaton condensate reaches its maximum value before oscillating back towards small 
values of 0. We denote hy t — —to the time when the inflaton condensate starts to roll away from its initial value 
(j) = 4>0i by f = —tm the moment when it reaches the value (f) = 4>m and by t = —ti the moment when it reaches the 
value (j) — see ([55)1 . Thus we have < ti < < to, see Fig. [5] Since here we neglect the expansion of the universe, 
the time evolution is symmetric around t = 0: as the inflaton rolls back towards small values of (f> after t — 0, it takes 
the same values (j) = 4>i, 4>m and (j)o at the times t = ti, tm and to respectively. For the modes with k > fcmax, the 
frequency squared in (|45)) is always positive. For the modes with k < /cmax, it is negative during a fc-dependent 

interval of time in the tachyonic regime. We will denote this time interval by when t < and [t^,^^] 



Note in this respect that a decrease of the modes in time can be viewed as an increase of the wave-functions in space, as the orientation 
of the coordinate t corresponding to the time-evolution of the modes is "opposite" to the one corresponding to the spatial profile of the 
wave-functions. Indeed, for the time evolution of the modes, the initial conditions correspond to only positive-frequency waves {a^ = 1 
and /3fe = in Il49ll i in the vacuum, while in general both positive- and negative-frequency waves are excited after one condensate 
oscillation when particles are produced. In the opposite way, in the quantum mechanical problem of tunneling through a potential 
barrier, both incident and reflected waves are present before the barrier, while only transmitted waves are present after the barrier. 
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when t > 0, see Fig. [5] Thus 



V"{^t)^0 (48) 




where (/)^. = 0(±tj,). The modes with k < kg, Eq. (|46p . have a negative frequency squared from the begining of 
preheating, so that = to for these modes. 

A convenient way of solving the mode equation, that is often useful in the context of preheating (see e.g. [lol.[26l|). is 
to use the WKB approximation when the frequency varies adiabatically with time, i.e. when \ujk\ <C and \u}k\ ^ 
(including imaginary values of ujk when < 0). Here however, we find that these conditions are satisfied only in the 
intervals of time [—tQ,—tm] and [tm,to\- Let us first consider the evolution of the modes in these intervals of time. 
For the modes with ko < k < k^i^^, see Eqs. (j46l47l) . we have t^^.(i) > for —to < t < — and for <t < to- The 
solution of Eq. (1451) can then be written as a superposition of positive and negative frequency waves 

Vkit) ^ -^^=exp f-z f dt'u^kit')) + -it=exp f dt' u^kit')) (49) 
for —to <t< and 

Vkit) ^^^ exp l-i f dt'iUkit')] + ^A^ exp ( i / dt' Wk{t') ] (50) 

yy2ujk[t) \ Jt+ J \/2uJk[t) 

for < t < to- The Bogolyubov coefficients at, f3k, (Xk and f3k are constant when the frequency evolves adiabatically 
with time. For the modes with k < fcmax, when the frequency squared is negative the solution of (j45p can be written 
as a superposition of exponentially increasing and decreasing solutions 

"'<''-™°--'-'(-/V''"'<''') 

for —t^ < t < —tjn and 

for tm ^ t < , where fi^ — ~uJk is positive. Again, the coefficients a^, bk, Ofe and bk are constant in the adiabatic 
regime. The second term in Eqs. (I51l52p is exponentially growing with time and corresponds to the tachyonic 
amplification of the modes. However, if the coefficient of the first, exponentially decreasing term is much larger than 
the coefficient of the exponentially increasing one, then the mode will start to decrease exponentially with time and 
this will explain the behaviour observed in Figs. |3] and ID 

The WKB approximation is not valid in the vicinity of the turning points t — ±i^, where the frequency vanishes, 
but there the oscillatory and exponentially evolving solutions can be matched, as described below (see also j2^), and 
this allows to follow the evolution of the modes through these moments of time too. However, the WKB approximation 
is not applicable around the entire crater of the Schrodinger potential, i.e. when —tm < t < tm, because the potential 
varies more abruptly there. Thus, in that region, it does not seem possible to find an approximate solution for the 
exact Schrodinger potential. In order to follow the evolution of the modes in that region, our strategy will be instead 
to find the exact solution for an approximate Schrodinger potential. To sum up, we will solve the mode equation (j45p 
for the modes with k < fcmax as follows: (i) use the WKB approximations (|49l - [5^ in the intervals of time [—to, —tm] 
and [tm,^o], (ii) match the solutions p9)) with (|5T|) and ((50)) with (l52t around the turning points t = for the 
modes with k > ko, (iii) find an approximate Schrodinger potential in the interval [—tm, tm] for which the Schrodinger 
equation can be solved analytically and (iv) match this exact solution with the WKB approximations ([STl [52l) at 

t — ^tm- 

We can apply the same procedure for the modes with k > /cmax, except that in this case we will use the oscillatory 
solutions ()49l I50|) in the WKB approximation in the entire intervals [—to,—tm] and [imjioli since the frequency 
squared is always positive for these modes. These oscillatory solutions will then be matched to the exact solution of 
the approximate Schrodinger equation in the interval [—tm,tm]- Note that, although the modes with k > fcmax are 
not amplified by the tachyonic instability, they can still be amplified by the non-adiabatic evolution of the frequency 
in the interval [—tm,tm]- We will see however that this effect is in fact negligible. 

The method outlined above allows one to follow the evolution of the modes and requires that we can indeed use the 
WKB approximation at t = ±tm, i-e. that the frequency varies adiabatically at these moments in time. This is very 
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well satisified for the modes with k <C fcmax, which dominate the second stage of preheating, because Cl^ oc V"'{4>) (j) 
vanishes a.t t = since T^"'(0m) — 0. Similarly, the WKB approximation is very good a.t t — for the 

modes with k ^ fcmax- However, it is not applicable for the modes with k w fcmax, because the frequency vanishes 
around t — ±tm for these modes. In principle, we could describe the evolution of these modes by finding another 
approximation of the Schrodinger potential that is accurate also outside the interval [—tm,tm]- However, this is not 
necessary, since we will see that our results are accurate even for k only slightly less or slightly bigger than fcmax, and 
the behaviour of the modes with fc w fcmax can be easily understood as an interpolation between these two regimes. 

It remains to find a Schrodinger potential which is a good approximation of —V" {<t>{t)) in the interval [— im,^m] 
and for which exact analytical solutions of the Schrodinger equation can be found. As discussed in the Appendix, a 
good choice is given by 

- V" m)) ^ U{t) = fcLx - (53) 



cosh [fit) 



with 



The above relation between the parameters does not follow from the only requirement that U{t) gives a good ap- 
proximation to —V". Indeed, we show in the Appendix that ^ should satisfy (j54l) in order to reproduce the known 
evolution of the zero mode Vk=o (x (p. We also show there that for the model (|14p with p = 2, the exact Schrodinger 
potential —V" appearing in the mode equation is given by (j531 154p with other numerical values of the parameters. 
As discussed in Section |TT1 this model does not lead to small field inflation, but in that case both the background 
evolution of the inflaton condensate and the wave-equation for the perturbations can be exactly solved analytically. 
For the model dST]), cj) is maximum with = at i = 0, so Eq. gS]) gives = 2^/* and therefore -V" = 11^2. This 
is reproduced by ((55)) for 

C^o = fcmax + 11V2~ 16.3 (55) 

where we have used (|47l) . Because ^tm >> 1, U{t) reproduces also the correct value fcmax of at = 0,„. The 
potential U {t) with the above values of the parameters is shown in Fig. [51 We see that it agrees very well with the 
numerical solution for —V" (0(t)) in the interval [— tm,im]- 
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FIG. 6: Comparison between the exact function ~V"{t) (in blue) calculated numerically and the analytical approximation U{t) 
(in red) defined in (|53|) . The agreement is very good in the interval [—t,n,t,n]. 

In the Appendix, we solve the Schrodinger equation for the potential U{t) and match the solution to the WKB 
approximations (1491152^ . In the next sub-section, we use these results to study the second stage of preheating for the 
modes with fc < fcmax- In sub-section IIV B|, we do the same for the modes with fc > fcmax- We then study the effect of 
the expansion of the universe on the second stage of preheating in sub-section IIV CI Finally, we generalize our results 
to other models of small field inflation in sub-section IIV Dl 









U(t) 
-v ' (t) 







A. "Tachyonic Oscillations" for k < fcmax 

Consider first the modes with fc < fcmax, still without the expansion of the universe. In this case, we have to match 
the exact solution of the approximate Schrodinger equation in the interval [—tm,tm] to the WKB approximations (|5ip 
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and (|52p. As shown in the Appendix, this leads to the transfer matrix 

ak\ X.- f ^ ~1 \ f "fc 



between the solutions (|5ip and (I52p in the limit fc ^ fcmax- We will see however that this approximation is valid even 
when k is only slightly smaller than fcmax- Here E ~ 2.15 is a numerical constant and 



r-*- J\V"U)\-k'^ 

Xk^ I Uk{t)dt^ I ^ ^ dc^ (57) 



is the "exponential phase" accumulated during half a period of the condensate oscillation when < and the WKB 
approximation is applicable. The lower limit of the integral above depends on the wave-number. The modes with 
fc < fco ~ ^/i(j)o - see Eq. PS|) - have a negative frequency squared from the beginning, so that (j)^ — (po for these 
modes. On the other hand, for the modes with k > ko, (jfj^ is the turning point where the frequency vanishes. For 
fc < fcmax, we have 0;!^ < 1 and Eqs. (|44l|48l) then give 0^ ~ fc/\/3. Thus 



-f= for k>^/i(j)Q 

^ (58) 
cj)(i for fc < V 3 (/)o 



for the modes with fc ^ fcmax- 

It is interesting to understand the behaviour of different modes when the condensate climbs back the potential 
from 0^1 towards (j) = (f)Q. For the modes with fc/fcmax <, e""^'', which includes in particular the zero- mode fc = 0, 
Eq. ([55)1 gives d^ ^ —e^^ bk and bk ^ e^^*" ak- Plugging into ([5^ . we then have Vkit^) ^ [bk — ak)/^/'2.Vlk which, 
from (j5ip . is of the order of Wfe(— t^). Thus we recover the fact that the low- momentum modes oscillate without any 
net growth during one oscillation of the condensate, because Vk oc </> for these modes. Indeed, for the modes with 
fc/fcmax <y e~-^'' , we have \dk/bk\ ^ e'^^'' and the solution ([52|) decreases exponentially during the whole interval of 
time [t,n, ] - Now we consider the modes with e~^'' <C k/k-a^sx ^ 1, for which the factor in e~^'^'= in ([56|) is neglibile. 
For these modes, gives \dk/bk\ ~ fc^/fcmax' which is still very large so that the solution ([5^ will again start to 
decrease exponentially with time, but not enough for this solution to decrease all the way until t = . 

In order to follow the evolution of the modes with fc > fco during a full oscillation of the inflaton condensate, it 
remains to match the eponential solutions (|51l52p with the oscillatory solutions (|49I50[) around the turning points 
t = ±t^, where the frequency vanishes and the WKB approximation breaks down. If we formally extend the domain 
of definition of the functions Vk{t), uJk{t) and ^kit) to the complex plane t, then the matching can be done by going 
around the real turning points along a contour in the complex plane which is sufficiently far away from the turning 
points that the WKB approximation remains valid along the contour [36j . This procedure leads to the transfer matrix 
(see [26| for more details) 



bk) \ e-^«^ ze^e. ) 



between the solutions (l49l) and (ISD) around t = — tm, where 



(59) 



Qk= / u:k{t)dt^ / ^ ^ dc^ (60) 

Jti J<t>o m 

is the " oscillatory phase" accumulated during half a period of the condensate oscillation when > and the WKB 
approximation is applicable. Proceding similarly around t = gives 

The factors in e"^'' in this transfer matrix are again negligible for k/k^^i^x > 6^^*= and will be set to zero in the 
following. 

Combining (15^1) and (pT|) . we get the full tranfer matrix between the solutions (|49p and ([501) during a complete 
oscillation of the inflaton condensate. We will denote by and (31 the coefhcients in after the j"^ complete 
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oscillation of the condensate (i.e. at the beginning of the (j + 1)*'^ oscillation) and by a^, ^ and f3l ^ the coefficients in 
(l50l) during the j^^ oscillation of the condensate. We then have the further relations aj, = (x'jT^ e"*®* and = e^^^ . 
We can then relate the Bogolyubov coefficients after and before the j**^ oscillation of the infiaton condensate as 

Because here we neglect the expansion of the universe, all the coefficients in this transfer matrix are the same during 
each oscillation of the condensate. The occupation number of the infiaton fluctuations rv\. = |/3^P after j complete 
oscillation of the infiaton condensate is then obtained by taking the j*'^ power of the tranfer matrix. With the initial 
conditions = 1 and — Q corresponding to vacuum fiuctuations, we finally get 
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E e^-^^ 

fc2 



2j 

[2 cos(2efe)]'(-'~^^ . (63) 



It is interesting to compare this simple analytical expression to the one obtained in |26| for the process of "tachyonic 
resonance" that takes place during preheating after chaotic infiation with trilinear interactions. In that case, the 
frequency squared of the modes becomes also periodically negative during some fc-dependent intervals of time, but 
the modes are only growing exponentially during these tachyonic episodes. In Eq. (j63p . 2Xk and 20^ correspond 
to the total phase of $7^ and ujk accumulated in the adiabatic regime during a complete oscillation of the infiaton 
condensate, see (|57p and (|60|). The expression (j63| is then formally similar to the one obtained in 26], except for the 
crucial pre-factor in /k"^^^ that takes into account the effect of the temporary exponential decrease of the modes 
that occurs in preheating after small field infiation. 

As in [2^, after more than one oscillation of the condensate, the modes with cos(20;i.) = are not amplified because 
of destructive interferences between the successive tachyonic instabilities. These modes form narrow stability bands, 
separated by wide instability bands where the growth of the modes is governed by the factor in k^ e^-^'' in Eq. (|63p . 
To calculate Xfe, we can use Eqs. (|43l44p in the second equality of (|57)) . For k <^ fcmax, the integral is dominated by 
its lower limit 0^ <C 1. This gives 



where we neglected higher powers oi cf) inside the integral. The value of 0^ is given in Eq. (|58p . For fc ^ fep, we 
can neglect the term in 0q in (|64l) . This gives 



e^^" « ( -— ) (65) 



A,, 



2 



max 



-2(^6-1) 



for ko <^ k <^ fcmax- Thus, for these modes, the spectrum of the occupation number (1631) after j oscillations of the 
condensate decreases with fc as fc to the power —2j{\/6 — 1). On the other hand, for k <^ fco, (pt — 4>o and Xk is 
approximately constant, so that n^, grows as fc*-' . In particular, the peak of the spectrum is located in between these 
two different regimes, at A: ~ fco- Note also form ([65t that fc/ fcmax ^ e~'^'' for fc ^ fco. This implies that only the 
modes with fc <C fco (i-e. well below the peak of the spectrum) have non-negligible factors in e"^"^*" in the transfer 
matrices (1551) and (pT|) . Therefore, as discussed in the paragraph below (|55|) . only the modes with fc <C fco decrease 
exponentially during the whole interval of time [im,to]- We will see that the expansion of the universe will modify 
this result. 

The spectrum of the occupation number of the infiaton fiuctuations after j oscillations of the condensate is shown 
in Fig. [71 where we compare the result obtained by solving numerically the mode equation to the analytical result 
(|63p in the range fco < fc. The aggreement is very good from the peak at fc ~ fco up to fc ^ fcmax- The left pannel 
shows the spectrum after j = 1 oscillation of the condensate, while the middle and right pannels are for j = 3. In the 
second case, we have narrow "stability bands" where cos(20fc) = in (|55)) . When the expansion of the universe is 
taken into account, the physical momenta k/a "move" inside the resonance pattern. This makes the band structure 
less distinct, but still present as we will see in sub-section IIV CI 



B. Non- Adiabatic Production for fc > fc, 



For the modes with fc > fcmax, the frequency squared remains always positive but the infiaton fiuctuations are still 
amplified in the interval of time [—tm, tm] where the frequency varies non-adiabatically with time. To study this case. 
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FIG. 7; Spectrum of the occupation number of the inflaton fluctuations (without expansion of the universe) after j complete 
oscillations of the condensate, for j — 1 (left pannel) and j — 3 (middle and right pannel). Here we took 00 ~ 10^^ so that 
fco — 0.017, see (146 |l . while fcmax — 0.87, see (|47|l . We compare the analytical expression (|63|l (in red) with the spectrum 
obtained by solving the mode equation numerically (in blue). The right pannel is a zoom of the middle pannel on a shorter 
range of momenta. 



we match the WKB solutions (l49l) and ([SO]) to the exact solution of the approximate Schrodinger equation in the 
interval [—tm, ^m]- As shown in the Appendix, we can then relate the occupation numbers of the inflaton fluctuations 



and 



after and before the j oscillation of the condensate as 



li = \D,\' + (1 + 2 \D,\^) + 2\Du\ V^+W ^ ni'^ (nf' 



l)cos0i: 



(66) 



where 0{ is a phase defined in the Appendix and 



\D,\ = 



sin (7rfcmax/M) 



sinh 



(67) 



Eq. (|66l) is formally identical to the one obtained in 10] for the non-adiabatic production of particles after chaotic 



inflation, but in that case \Dk\ 



^fc /fc. where fc* is the typical momentum amplified by the process. 
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FIG. 8: Spectrum of the occupation number of the infla- 
ton fluctuations (without expansion of the universe) after 
3 oscillations of the condensate, for k ~ fcmax — 0.87. 



FIG. 9: Evolution with time of the (comoving) occupa- 
tion number nu for modes with k ~ fcmax, for the model 
(|37p with expansion of the universe and v — 10^^ GeV. 
The evolution is shown during the first 6 oscillations of 
the condensate, where the scale factor varies from a ~ 1.5 
to a ~ 1.7 (with a = 1 at the beginning of preheat- 
ing). Thus the maximum comoving momentum ampli- 
fied by the tachyonic effect, fcmax — 0.87 a, varies from 
fcmax — 1.3 to fcmax — 1.5. We show the evolution of the 
modes with comoving momenta (from top to bottom): 
k = 1.2 (blue), k = 1.5 (red) and k = 1.7 (black). 



(X e '^^/^^ in 



For k > fcmax, \Dk. 

as k decreases towards fcmax, 



66|) and the non-adiabatic production is very inefficient. On the other hand, 
Dk\ in (j67|) can become smaller than one and the efficiency of the non-adiabatic 
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production is enhanced. This is to be contrasted with the non-adiabatic production after chaotic inflation studied 
in [l3|, where |-Dfcp = e"'^''^/'^* < 1. However, as discussed in the paragraph before Eq. ([55)) . the matching method 
and therefore Eq. (j66l) are not valid for k w fcmax- In this limit, there is a smooth transition between the tachyonic 
amplification for k < fcmax and the non-adiabatic production for k > fcmax. Indeed, when k increases towards fcmax, 
the frequency squared becomes negative in shorter and shorter intervals of time around itm and the efficiency of 
the tachyonic amplification decreases considerably, whereas the efficiency of the non-adiabatic production increases 
when k decreases towards fcmax- This is illustrated in Fig. [8l where we show the spectrum of the occupation number 
of the inflaton fluctuations for k « fcmax after 3 oscillations of the condensate. Note again the stability bands where 
destructive interferences occur. Their positions depend on 0^. in (j63l) for k < fcmax and on Oj. in (|66p for k > femax- 
Note also that, even if the non-adiabatic production is enhanced for k « fcmax, it is still much less efficient than the 
tachyonic amplification of the modes with fc <C fcmax discussed in the previous sub-section. 

For the modes with fc ^ fcmax, particle production occurs on the physical (un-rescaled) time scale tm ~ ™min ^ 
during which the expansion of the universe is negligible. For these modes, the main effect of the expansion of the 
universe is that the physical momentum k/a varies from one event of particle to the other. In that case, as in [lo| . 
the physical momenta move inside the resonance pattern of Fig. [8] during successive oscillations of the condensate. 
This implies in particular that the occupation number of a given comoving mode varies in a different way at different 
events of particle production. This is illustrated in Fig. ^ where we show the evolution with time of the (comoving) 
occupation number in an expanding universe for modes with fc « fcmax- Because of the last term in (|66p. the occupation 
numbers can decrease during one oscillation of the condensate, see e.g. the mode with fc = 1.5 in Fig.[Sl Furthermore, 
since k/a decreases with time whereas k-^aax/a = \/ —V" {(t)m) remains constant, more and more modes enter the 
regime with k/a \/ ~y" {(f'm) where the tachyonic amplification becomes more and more efficient. This occurs for 
instance for the mode with fc = 1.2 in Fig. [9l which start with fc/a « y^—V"{(j)m) and ends with k/a <^ \/—V"{(j)m)- 
In the following sub-section, we discuss the effects of the expansion of the universe for the modes with fc <C fcmax, 
which dominate the second stage of preheating. 



C. Tachyonic Oscillations in an Expanding Universe 

We now discuss how the process of "tachyonic oscillations" studied in sub-section llV Al for the modes with fc ^ fcmax 
is modified by the expansion of the universe. For these modes, the main effect of the expansion of the universe comes 
from the fact that the amplitude of the condensate oscillations decrease with time due to Hubble friction. The effect 
of Hubble friction on the evolution of the condensate was discussed in Section |lTl where we obtained the estimate ((22|) 
for the minimum value 4>j of the condensate when it climbs back the potential after j > 1 complete oscillation(s). For 
the model (|57|. it reads 

\ 1/4 



^jyi'^ij — • (68) 

Here we continue to use the rescaled variables (|42p. but now we write the "tildes" explicitly to avoid confusion. In 
these variables, the initial value (|38p of the condensate at the beginning of preheating reads 

4>j=0 =<t>e^ —f^ • (69) 

V24 TOpi 



Thus is of the same order of magnitude for different values of j > 1, but (^j>i ^e- 

First of all, as discussed in [2^, Hubble friction may be so important that the condensate never goes back in the 
tachyonic region </> < after his first rolling down. Using Eqs. ([55)1 and i^i > <^rn for 

— > 10^1 . (70) 

TOpi ^ 

This is illustrated in Fig. I10[ where we show the evolution with time of the inflaton condensate during its first 
oscillations for different values of v/m-pi- When the condition (170)) is satisfied, the process of tachyonic oscillation 
does not take place. The tachyonic amplification of inflaton fluctuations during the flrst stage of preheating is also 
rather inefficient for such high values of w/mpi. In that case, the main decay channel for the inflaton condensate is 
the decay into other fields, either perturbatively or non-perturbatively. The decay of the inflaton into other fields will 
be discussed in the next Section. 

We also saw above that, for the model ([57]) . preheating ends after the first stage (i.e. in less than one oscillation of 
the condensate) for v/mpi <^ 10~^. Thus we now discuss the effects of the expansion of the universe on the second 
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stage of preheating for 10 ^ < v/mpi < 10 ^. 
during the j"^ oscillation of the condensate 



In order to do so, it is convenient to define the analogue of (|57l [58]) 



where 



XI = 



0" Vl^"l - '^Va^ 



\/3 a 



for 
for 



(71) 



(72) 



The scale factor a does not vary much during the first few oscillations of the condensate after the first one. We will 
take it to be constant and equal, say to its value after the first complete oscillation of the condensate. Note also that 
we consider the modes with k ^ fc,nax — 0.87, so that k/a ^ 1. Proceding as in (|64l [65]) . we then have 



constant 



for ^ > 

for I < 73^, 



(73) 



where we work at the level of an order of magnitude estimate 
depends on time. 

To see this, consider the j*^ oscillation of the condensate, where 



The main effect comes from the fact that Xf, now 



' starts from 



oscillates around the minimum 



of the potential and climbs back the potential up to We can then use the transfer matrix (|56|) with Xk replaced 
by X^i,~^, fc^ replaced by the physical momentum k^ /a? and E/kf-^^^ ~ 1. On the other hand, in the transfer matrix 
([5T|). Xk should be replaced by Xi. This leads, instead of to the estimate 




-I 2 



-.XI 



4j 



2X° 2X1 



AXl 



(74) 



for the occupation number of the inflaton fluctuations after j complete oscillations of the condensate. Indeed, during 
the j*^ oscillation of the condensate, is first amplified by a factor of e^-^k when the condensate rolls from 4> = 4>j~i 
to the minimum of the potential, and then by a factor of e^-^k k'^/a^ when the condensate rolls back up to (/) = 4>j. 
As before, the factor in k'^/a'^ <C 1 results from the fact that the modes start to decrease exponentially when the 
condensate rolls back towards small values of (j). 

Note that, as in sub-section IIV Al Eq. ([71]) holds only for the modes which eventually increase exponentially when 
the condensate rolls back towards small values of <j> (instead of decreasing exponentially during this whole stage), 
which dominate the second stage of preheating. These are the modes for which the factors in e~^'' in the transfer 
matrices (1551) and (pT|) are negligible, i.e. the modes with k/a > e~^i'. Using (|73p with k/a <C 1, this condition 
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reads k/a > ■ Since 0j in is of the same order of magnitude for different values of j > 1, Eq. ([7^ holds 

for the modes with k/a > (j)^ . On the other hand, Eq. (|74p is not valid for the modes with k/a < 4''f^, which 
decrease exponentially during all the time that the inflaton rolls back towards small values of (f>. These modes behave 
as the zero-mode Vk oc (j) and oscillate without any net growth during the second stage of preheating. One can 
check that these include in particular the modes of the Hubble scale at the end of inflation, i.e. the modes with 

k/a ~ Hg/^imin ~ ^0 ^ 4>'i^ in rescaled variables . Indeed, we already saw in Section Hill and Fig. [3] that these 
modes follow the condensate velocity, Vk oc (p. 

In the following, instead of the occupation number, it will be more convenient to consider the variance per logarithmic 
momentum interval, fc^ \vk\'^- Note from (|4T|) that rik ^ |wfcP k/a. Then, using (|74)) and k — 2fc/mmin, we have 



k 



■ n e'^^ (75) 



after j complete oscillations of the condensate. Because of the two different regimes ((75)) for each factor in e"^^- in (|75p , 
the spectrum has a different shape in each interval -v/S^i-i < k/a < V3(f)i. For instance, for < k/a < \/3^i, all 
the factors in 6'^'= are constant, except (since we saw above that i^o ^ "^i^)- In that case, the spectrum increases 
as k to the power Aj + 2 — 2\/6. On the other hand, for k/a > \/3(f>j, all the factors in e'^'-- in ([75)1 decrease as k~^. 
In that case, the spectrum decreases as k to the power — 4j(\/6 — 1) + 2. In particular, the peak of the spectrum is 
located in between these two regimes, at A:*/a ~ V^'Pi with 1 < i < j. Using (I68p . this gives k^/a ~ {v/mpi)^^^ . In 
terms of the un- rescaled momenta p2)) . we have 



a \mpi 



mmin (76) 



for the characteristic physical momentum amplified during the second stage of preheating. Note that this is much 
larger than the characteristic momentum amplified during the first stage of preheating, fc*/a ~ H^., so the spectrum 
is shifted towards the UV by a very large amount during the second stage. 

From l|75p. we can estimate the variance of the infiaton fiuctuations tpk — a~^/'^Vk after j oscillations of the 
condensate as 

e-^e-^. n^^"- (77) 

where the RHS can be evaluated at the position of the peak, k^/a w (v/mpi)^^'^. Using (j39p and (|73l) . this gives 

MW V \-(v^-i)J+i/2 



mpi 



(78) 



We can now estimate the number of oscillations of the condensate after which (ScjP) ~ and preheating ends. In 
([78]) . M depends on v through the normalization of the CMB perturbations, Eq. ([rU]) with p = 4 for the model (|37p . 
For the range 10^^ v/mpi <^ 10~^ that we are interested in, the number of efolds before the end of inflation when 
cosmological scales leave the Hubble radius is of the order A^* w 45 - 50. We then find that preheating in the model 
(|37ll ends after only 



10 



logio (mpi/v) 



(79) 



oscillations of the condensate. This varies from j ~ 2 for v/nip\ — 10 ^ to j — 5 for v/nip\ = 10 ^. Thus, for 
v/mpi <^ 10^^, preheating is still very efficient despite the temporary exponential decrease of the modes. 



This is in contrast to the case without expansion of the universe, where we saw that only the modes with fc <C fco ~ </'0 decrease 
exponentially during all the time that the condensate rolls back towards small values of <j)- This is because in the absence of Hubble 
friction, the condensate rolls back to its initial value 0o ^ 4'1 after each oscillation. The second (exponentially growing) term in 1521 1 
has then much more time to catch up with the first (exponentially decreasing) one, so that many more modes eventually increase 
exponentially when the condensate is rolling back. 
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FIG. 11: Spectrum of inflaton fluctuations fc^ IV'feP a-t the end of the hnear stage of preheating (when iScf)^) « u^) in the model 
(f37)l for ij/mpi = 10"^ (top left), 10"^ (top right), 10"" (bottom left) and 10"^ (bottom right). Preheating ends after j ^ 6, 4, 
3 and 2 oscillations of the inflaton condensate, respectively. The spectra are normalized in such a way that the peak amplitude 
is of order one at the end of preheating and they are shown as a function of the physical momentum k/a in units of mniin/2. 



The spectrum of the inflaton fluctuations calculated numerically at the end of the second stage of preheating is 
shown in Fig. [11] for different values of v/mpi. The results are in very good agreement with the analytical estimates 
(j76p and (j79p . At momenta larger than the peak, some band structure remains due to destructive interferences, 
despite the expansion of the universe. 



D. Generalization to other models of small field inflation 



The method that we developed above to study the second stage of preheating in the model ([57)) is generic and can 
be applied to other models of small field inflation. To illustrate this, we now generalize our results to the models ([H]) 

V = - J) (80) 

with any p > 2. As discussed in Section [TTJ p > 2 is the condition to have a small field inflation model. In this case, 
the inflaton mass at the minimum of the potential is given by 

M2 

rUn^in = VV"{v) = V2p . (81) 

V 



The maximum of —V"{(p) is located at 

. / p-2 



"^'^ \2{2p-l) 
and the maximum momentum amplified by the tachyonic effect is 

h 



V (82) 



(83) 
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As before, Hubble friction makes the second stage of preheating inefficient if <^ (pi , where the amphtude of the 
condensate 4>i after one complete oscillation is now given by ([25)1 with j — 1 and X — 2p for the model (|80p . Comparing 
with (|82|) . we find that the condition (f70| for the tachyonic effect to be inefficient remains approximately the same 
for any p. Similarly, the condition p6p for preheating to end during the first stage gives v/mpi <^ 10^^ — 10~^ for any 
value of p. Thus, as before, the second stage of preheating in the model (IMl) take place for 10~^ <; v/mpi ^< 10^^. 

The second stage of preheating is again dominated by the regime of tachyonic oscillations for the modes with k <C 
kmax- This can be studied analytically for any p by a straightforward generalization of what we did in sub-sections lIV A1 
and IIV Cl The efi'ective Schrodinger potential —V"{(j)) in the mode equation (|45p can again be approximated by the 
RHS of ([53]) in the region ip > (pm- This is illustrated in Fig. [12] for p = i and p = 8. The parameters fcmax, Uq 
and jJL take different numerical values for different values of p, but they are still related by Eq. ([54]) . This implies 
in particular that the WKB approximations ([?T|) and (|52p are still related by the transfer matrix ([55]) . with the 
characteristic suppression in fc^/fc^^ax ^ 1 due to the temporary exponential decrease of the modes when </> < 0. 




FIG. 12: Same as Fig. H but now for the model (gnj with p = 3 (left) and p = 8 (right). 
On the other hand, the rate of tachyonic growth depends on the value of p. Defining 

_ v/2p(p-l) 



(84) 



and proceding as in (j73p . we find the estimate 



for 
for 



a ^ 



(p-2)/2 
3 

^(p-2)/2 



for k/a 1, where the "tildes" denote the rescaled variables p2|) and 



V 

mpi 



(85) 



(86) 



is the amplitude of the condensate after j > 1 complete oscillations. As p increases, for a fixed ratio of v/mpi, 
the effect of Hubble friction on the amplitude of the condensate (|86|) increases, so that the tachyonic growth of the 
low-momentum modes occurs on a shorter range of field values. However, the rate of growth of these modes increases 

also with p and their net tachyonic growth e^i ^ {mpi / v)^''^^^^^^ increases slightly as p increases. 

We can then proceed as we did below Eq. ((75|) . At the end of preheating, the peak of the spectrum of the infiaton 



fluctuations is now given by fc*/a « (p^ 



(p-2)/2 



In terms of the physical variables, this gives 



p-2 
2p 



(87) 



-(xp-l)2j + l] 



Therefore, in terms of m^in, the peak of the spectrum moves towards the IR as p increases. This is because, as we 
noted above, the rate of growth of the low-momentum modes increases when p increases. The variance of the infiaton 
fiuctuations grows as 

p-2 r 

\ mpi 

after j complete oscillations of the condensate. As p increases, the rate of growth of the fiuctuations increases, but 
the ratio M'^/v'^ fixed by the CMB normalization (ITUl) decreases. For any reasonable value of p (say 3 < p < 20), we 
then find again that preheating ends in between 2 and 6 oscillations of the condensate for 10~^ <^ v/mp\ <; 10~^. 
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V. OTHER NON-PERTURBATIVE DECAY CHANNELS, PERTURB ATIVE DECAY AND REHEAT 

TEMPERATURE 

We now study the decay of the inflaton into other fields, for a generic model ([l] [2]) of small field inflation. We first 
consider other non-perturbative decay channels for the condensate and show that they are usually negligible with 
respect to the decay channel into inflaton fluctuations that we studied in the previous sections for v/mp\ <^ 10^^. 
Preheating is then dominated by the mechanism that we studied above, leading to the very quick decay of the inflaton 
condensate into large and non-thermal fluctuations of itself. As usual, this is only the first stage of reheating. The 
inflaton fluctuations must then decay into other degrees of freedom and the universe must eventually thermalize. 
When the minimum of the inflaton potential 4> — v corresponds to a fixed point of internal symmetries, the inflaton 
fluctuations may decay into massless fields and thermalization may be relatively fast. On the other hand, when cf) = v 
is not a fixed point of internal symmetries, the non-zero VEV of the inflaton at the minimum of the potential provides 
a large effective mass to the fields coupled to it. As a result, the inflaton can decay only into the fields to which it 
is weakly coupled, which therefore delays the decay. We will thus consider these two different cases separately in the 
following two sub-sections. In each case, we determine the reheat temperature and the resulting number of efolds iV, 
before the end of inflation when cosmological perturbations leave the Hubble radius. We then determine, for different 
values of p, the explicit relation (|TOl) between M and v that follows from the normalization of the CMB anisotropics, 
and the region in the parameter space for which preheating ends in less than one oscillation of the condensate. 



A. Restored symmetries at the minimum 

Let us first consider the case where the minimum of the infiaton potential (p = v corresponds to a fixed point of 
internal symmetries. This is not a very frequent situation in small field models, but it arises for instance in the model 
of 33]. In that case, the inflaton can couple to massless flelds at (j) — v. We can consider for example the interaction 
term 

^{cj^-vfx' (89) 

of the inflaton with another scalar held %, which is massless when (j) = v. 

In that case, if the inflaton condensate did not decay during the flrst stage of preheating - i.e. v/mp\ J> 10^^ — 10^^, 
see Section [nil - it can lead to the non-perturbative production of x-particles when it oscillates around the minimum 
(j) = V. If the ^-particles are relatively long-lived, their number density can grow due to parametric resonance [loj . 
The condensate oscillates with an initial amplitude of the order of $ < d and a frequency of the order of l/mmm- 
The efficiency of the resonance is governed by the dimension-less parameter 



^2 
,2 



«rT7I (90) 



where we used mmin ^ A/^/w. The relation ()10[) shows that the resonance can be broad {q > 1) or narrow (q < I), 
depending on the value of and p. For p = 3 or 4, we have usually g ^ 1, which is the most efficient case. In that 
case, the occupation number of the v-particles grows as e^'^'=^ after j oscillations of the inflaton condensate, where the 
maximal value of /ife is /ifc ~ 0.28 [lO|. This growth is clearly negligible with respect to the mechanism of tachyonic 
oscillations that we studied in sub-section IIV Dl 

Another possible mechanism for preheating is called instant preheating 37]. This may occur if x-particles are 
produced non-adiabatically when (p = v and then decay into other particles in less than one oscillation of the inflaton 
condensate. Again, this requires that the condensate did not decay before reaching the minimum of the potential, 
v/mpi ^ 10^^ — 10^^. When it crosses the minimum, it produces massless x-particles with a number density [s^] 



(5 4>Sj 



3/2 

— (91) 

o vi - 

via the interaction ([89)) . Here ~ is the velocity of the condensate at the minimum (j) = v. After they have been 
produced, the effective mass of the x-particles increases as — g — v\. If these x-particles decay, for instance into 
fermions with an interaction term 



hxipi', 



(92) 
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then their decay rate is given by 

h'^ h'^g\(l}—v\ 
^x^^ i: = = ^ (93) 



which is maximum when |0 — ~ t; is maximum. For the mechanism of instant preheating to be efficient, this decay 
should occur in less than one oscillation of the inflaton condensate, i.e. T ^ rrirain- This condition gives 

/.^>^>H!^ (94) 
g\<p-v\ g 

where we have used mmin ^ M'^/v and |(/) — <^ v. The condition (IMl) may be easily satisfied for p = 3 or 4, see 
Eq. 

In that case, once the x-partides have been produced when (f> = v, they decay after a time At ^ l/F, when the 
value of the inflaton condensate is 4>d — v ~ At <py . This gives 

(^)" (-) 

when the x-particles decay. The decay produces V'-particles with a typical energy = where = g \(l)d — v\. 

The energy density of the V'-particles p^, — n.^ E.^ compared to the energy density of the inflaton condensate — 0^/2 
is then given by 

^ = !k:!k«0.05^ (96) 

where we have used Eq. (1911) . Thus this mechanism can be very efficient for suitable interactions between the different 
fields. For v/mp\ ^ 10^^ — 10^^ and h <^ g^, it can compete with the tachyonic decay into inflaton fluctuations if 
the minimum cf) — v is indeed a flxed point of internal symmetries. The last condition is necessary because otherwise 
the non-adiabatic production of x-particles is negligible from the beginning, see the next sub-section. A particularly 
interesting realisation of the instant preheating framework has been recently provided in [s^l where inflation occurs 
along the LLe D-flat direction in the MSSM. After the end of inflation, the inflaton oscillates around a point of 
enhanced gauge symmetry. Tachyonic preheating would end in a large (w 10) oscillations while in fact instant 
preheating implies that the inflaton loses enough energy in the flrst oscillation never to go back to the tachyonic 
region. In the subsequent oscillations, instant preheating transfers the energy from the inflaton to the SM quarks, 
leptons and their supersymmetric partners. 

Note also that, even in the more generic case where the dominant mechanism for preheating is the tachyonic decay 
into inflaton fluctuations, fluctuations of other flelds may then be quickly amplifled too via their interactions with the 
inflaton fluctuations ^39.] . Thus in general, at the end of preheating, the universe is dominated by large fluctuations 
of the inflaton and other flelds coupled to it. Since the inflaton fluctuations are massive, they become quickly non- 
relativistic and dominate the total energy density before they decay. This subsequent decay of the inflaton can be 
described by the perturbative theory of reheating pj|] . We now consider this stage of perturbative decay. 

When (f) = V is a flxed point of internal symmetries, the inflaton can couple to massless flelds and quickly decay into 
these. Trilinear interactions allow for the inflaton to decay completely. This is provided for instance by an interaction 
term 

h\(t)-v\^^ (97) 

with fermions ip, which are again massless for (p = v. The decay rate of the inflaton fluctuations S4> — (j) ~ v, whose 
effective mass is ^ « ?Timin, then reads 



When the inflaton fluctuations are large, (S(f)'^) ~ v'^, they contribute to their own effective mass which may then differ from ^yV"{v). 
This however does not change the order of magnitude. 
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In small field models, mmin ^ M^/w can be much larger than the Hubble rate during inflation, He ^ M^/mpi. If the 
coupling constant h is not too small, T > Hf. and the inflaton then decays in less than one Hubble time after the end 
of inflation. Assuming that the decay products quickly thermalize, the reheat temperature is then almost of the order 
of the energy scale during inflation 

Tfl - M . (99) 

We will call this case "instantaneous reheating". We will see that the situation can be very different when (p = v 
is not a flxed point of internal symmetries. Note also that Eq. ()99|) implies that inflation should occur at a rather 
low energy scale in order to satisfy the usual upper bound on Tr to avoid the thermal over-production of gravitinos, 
M ^ 10^ GeV or so. 

Let us now determine the explicit relation between M and v from the normalization of the CMB anisotropics ([TU]) 
in the case of instantaneous reheating. In general, the number of efolds before the end of inflation when observable 
scales leave the Hubble radius is given by 

where Minfl — 6.6 10^^ GeV, see ([9]). Here V^nfi = M"^ is the energy density at the end of inflation and pioh the energy 
density when the universe becomes radiation-dominated. When reheating is instantaneous, pich = Mnfi so that 

~ 63.9 + In f -r^ ) . (101) 

V "Mnti J 

Inserting this into pH)) gives the explicit relation between v and M. This is shown in Fig. [T3] for X — 2p and 
different values of p. The almost horizontal lines that are roughly superimposed at the top of the plot show the upper 
bound p6l) on v below which preheating occurs in less than one oscillation of the condensate, for the same values 
of p. For p — 5, preheating occurs in less than one oscillation of the condensate only for small energy scales during 
inflation, M <^ 10^ GeV. This upper bound increases with p, e.g. M <^ 10^ GeV for p = A, which is similar to the 
gravitino bound, and M <^ 10^^ GeV for p = 5, which corresponds to the intermediate scale of SUSY breaking in 
gravity-mediated scenarios. 
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FIG. 13: Relation between v and M from CMB normalization, for p = 3 (black), 4 (blue), 5 (green) and 10 (red), in the case 
of "instantaneous reheating". The almost horizontal lines that are roughly superimposed at the top of the plot show the upper 
bound (|36p on v below which preheating occurs in less than one oscillation of the condensate, for the same values of p. 
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B. Broken symmetries at the minimum 



We now consider the case where (f) — v is not a fixed point of internal symmetries. We will first consider the decay 
of the inflaton into fermions ijj and another scalar field x via the renormalizable interaction terms 

h^b^P^ , Y'^'x'- (102) 

These interactions arise naturally if = is a fixed point of internal symmetries, which may help to ensure the 
fiatness of the infiaton potential during inflation. At the minimum of the inflaton potential, the non-zero VEV v of 
the infiaton provides an effective mass to ip and x via the interactions (|102l) . = hv and = gv. 

As in the previous sub-section, let us first check if the interactions (|102l) can lead to non-perturbative decay channels 
for the inflaton condensate that could compete with the decay into inflaton fluctuations that we studied in Sections IIIII 
and IIVI We will consider the non-perturbative production of the scalar field x since the non-perturbative production 
of fermions would be similar but limited by Pauli blocking. Let us first suppose that any bare mass of the x-field is 
negligible compared to the contribution coming from the interaction (|102p . so that the effective mass of x is simply 
™x — 9 4>- Efficient production of x-particles may occur when the effective mass varies non-adiabatically with time. 



I.e. 



;> 1 . (103) 



From Eq. (|20p. we have p oc 0^/^ for p v. This implies that m^/ml is maximum at small field values for p < 4 
and at = u for p > 4. At the beginning of preheating, p = p^, Eq. (|103l) gives 



l [p-l + (p-2)iV.]-^^-^)/(^-^) A4fl 

9 ^/6(p-l)(''-^)/(^-i) "^li ^ ^ 

where we have used Eqs. ([5]), ^ and (fTO|) . This shows that ^ at small field values, except for very small 
values of the coupling constant g. However, for such small values of the contribution g(j) to the effective mass of x 
is smaller than the Hubble rate during inflation. In order for x to remain massive during inflation, we must include a 
bare mass > iJ, which in that case is not anymore negligible with respect to g (p. When this is done, we find that 
^ dX 4> — (pe for any value of g. Therefore, there is no non-adiabatic production of x-particles for p < 4. 

For p > 4, m^/m^ is instead maximum at the minimum of the potential, (p = v. We have p ~ y/2M^ there, so 
Eq. (fins)) gives 



ml 



— . (105) 

gv 



For 5^1, this is of the order of A/^/w^ which, from Eq. (1101) . can be larger than one only for p > 4 and v much 
smaller than mpi. However, as discussed in Section IHIl for such small values of v the inflaton condensate decays 
before it reaches the minimum of the potential, so that there is no non-adiabatic production of ^-particles at 4> = v 
in that case either. Therefore, the only case where such a non-adiabatic production is possible if for p > 4 and 5 <C 1, 
where may vary non-adiabatically around (/> = v if the inflaton condensate did not decay during the flrst stage 
of preheating. As in the previous sub-section, either parametric resonance or instant preheating are a priori possible 
in this special case. However, the condition ([M|) is not satisfied in that case, so the x-particles can only decay after 
several oscillations of the condensate and instant preheating does not take place. As for the production of x-particles 
by parametric resonance, it is again negligible with respect to the tachyonic ampliflcation of the inflaton fluctuations. 
We therefore conclude that, when (p = v is not a flxed point of internal symmetries, preheating is in general completely 
dominated by the tachyonic effect that we studied in Sections IHIl and ITVl 

As in the previous sub-section, we now turn to the perturbative decay of the inflaton after preheating. Since we 
consider inflaton fluctuations around the minimum, (p = v + Sep, the interaction term g^ (p>^ with the scalar fleld 
X now includes a trilinear coupling g'^vS<px- Thus both interactions in (|102l) allow for the complete decay of the 
inflaton. The corresponding decay rates are given by 

_ h'^m^ g^v"^ 

when the fermion tp and the scalar x are much lighter than the inflaton (p, and when H <C (which is always 
satisfied for H ^ T, when most of the decay occurs). The mass of the inflaton is again of the order of the curvature of 
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the bare potential at the minimum, « mmin ~ M"^ jv. On the other hand, as we already mentioned, when (j> = v 
is not a fixed point of internal symmetries, the interactions (|102|) provide an effective mass to ip and % at = 
TO^ = hv and = gv- However, the decay of the inflaton is kinematically possible only if its mass is larger than at 
least twice the mass of the decay products, mip/2 > m^^m^. This requires 



,9<^ (107) 
2 V 



so that Eqs. (jl06p give 



to5 



Thus the decay is possible only for sufiiciently small coupling constants, which in turn delays the decay. 

This kinematical blocking of the inflaton decay when (j) = v is not a fixed point of internal symmetries occurs 
because the non-zero VEV v of the inflaton at the minimum provides a large effective mass to the flelds coupled to it. 
On the other hand, if some fields couple to the inflaton through derivative interactions only, these fields may remain 
massless despite the non-zero VEV of the inflaton. One may therefore wonder if derivative interactions could provide 
the fastest decay channel for the inflaton. The inflaton may have such interactions for instance with a gauge fleld , 
with the trilinear couplings 

where F^i^ is the gauge field strength, F^'^ = e^'^P'^ Fpa/2 is the dual field strength and A is a mass scale such as the 
cutoff scale or the symmetry breaking scale v. Notice that the coupling involving the dual field strength implies that 
(j) should be a pseudo-scalar and therefore the inflaton potential should be an even function. The first and second 
couplings in (|109p are natural if the infiaton </) is a modulus or an axion, respectively. The corresponding decay rates 
are given by 

where, in the last inequality, we have taken the condition A > w. When the interactions (|109p arise as the leading term 
in a series expansion in powers of (j), this condition should be satisfied for the series to be in a controllable regime. For 
instance, if A is the cutoff scale of the effective theory, w < A is required for the model of small field infiation to lie 
within the regime of validity of the effective theory. Even if we had t; ^ A, the first (modular) coupling in (jl09l) with 
(j) = V -\- 5(t> includes a contribution v/{AK)F^^ F^'^ to the kinetic term of the gauge field. Rescaling to the canonical 
gauge field then leads to the same bound (|110l) for the decay rate of d4>. This does not arise for the second (axionic) 
coupling in (|109|) because F^i, F^'' is a total derivative. One could then imagine a situation where the inflaton has 
only such an axionic coupling and the decay rate would be greater than the bound in (IllOp . However, as we just 
argued, this is very unlikely. 

Interestingly, the upper bounds (|108[ IllOl) on the decay rates for each interaction are all of the same order of 
magnitude. In the following, we will assume the presence of at least one such interaction with a coupling constant 
that maximizes the decay rate. Using w Wmin ~ M'^/v (see e.g. ([5T|) for the model ([SO)) ), we then have 

r«— (111) 

for the decay rate of the inflaton. Compared to the Hubble rate at the end of infiation He = M^/ (\/3topi), this gives 

r Af4 mpi , , 

— « — ; . (112) 

Eq. ((To)) shows that F <C i?e in large regions of the parameter space. In that case, the universe enters a long 
matter-dominated stage after preheating, until the non-relativistic infiaton fluctuations decay completely. 

Let us first determine the conditions on the parameters for which, despite the suppression of the decay rate, the 
inflaton can nevertheless decay in less than one Hubble time after the end of inflation. This occurs if F > He- Using 
([T2l) in ([TT2]) . this gives 

A//"4 / M \ -2(3p-10)/p 
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where C is defined in (|TT|) . When this condition is satisfied (and assuming again that the infiaton decay products 
quickly thermahze), reheating is "instantaneous": ^ M and the number of efolds is given by (|10ip . As discussed 
below Eq. ([T2|) . we have usually C > 1 and M < Mi^R- This implies that the condition (|113|) cannot be satisfied for 
p < 10/3 ~ 3.33. Thus reheating is never instantaneous for p — 3. For larger values of p, the condition (|14p with 
(|111 llOip gives an upper bound on M below which reheating is instantaneous. This upper bound on M is shown in 
Fig- El for different values of p. For p — 4:, reheating is instantaneous only if Tr ^ M < 500 GeV, which is very small 
but can still be compatible with electroweak baryogenesis. For larger values of p, reheating is instantaneous for larger 
values of M, e.g. M ^ 10" GeV for p = 6. 

Let us now consider the other case, when T < He- In that case, reheating is delayed until the Hubble rate H drops 
below r and the reheat temperature is then given by Tr « 0.5 \/T mpj. Using (|llip . this gives 



A'P /mpi\i/2 



C ~ I -^T— I M 



'PI 



infl 



(114) 



where we have used in the second equality. The number of efolds is now given by (llOOp with prch — 3mp[ F^ 
Using (fTTn[T2)) . this gives 



7V, + ^lnC: 
op 



63.9 



10 , f M 
— In 

3p VMnfl 



(115) 



From this, we can calculate N^, and C as a function of M and then use (|114p to calculate Tr as a function of M. 
The result is shown in Fig. [15] for different values of p. The curves for different values oi p > 10/3 are superposed 
for sufficiently low M, because reheating is then instantaneous and Tr ~ 0.5 M independently of p. Then, when M 
becomes greater than a certain value that depends on p (the value shown in Fig. 1141) , reheating starts to be delayed 
and Tr increases with M more slowly than Tr (x M. In fact, Eq. (|114p shows that, for p > 5, the reheat temperature 
even decreases when the energy scale of inflation increases! This is because, in that case, the requirement that the 
perturbative decay of the inflaton is kinematically possible leads to a more stringent upper bound on F when M 
increases, so that reheating takes a longer time to complete. Meanwhile, for p = 3, we saw above that we have 
F <§C He so that the reheat temperature is typically much smaller than M. 

Finally, we show in Fig. [1^] the relation between v and M for different values of p. In each case, we only show 
the values for which v < mpi and Tji > GeV, and the dot indicates the values of M and v below which reheating 
is instantaneous. The almost horizontal lines that are roughly superimposed at the top of the plot show the upper 
bound ([55]) on v below which preheating occurs in less than one oscillation of the condensate, for different values of p. 
For p = 3, since the decay rate is very small when ip = v is not a fixed point of internal symmetries, preheating must 
occur in more than one oscillation for the universe to reheat before BBN. For the other values of p, all the cases are 
possible depending on the energy scale of infiation: preheating in less or more than one oscillation and instantaneous 
or delayed reheating. Each case will lead to different predictions for the cosmological consequences of preheating after 
small field infiation, such as the production of primordial black holes and gravitational waves. 
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FIG. 14; Maximum value of M (in GeV) for which the con- 
dition (|113p is satisfied, as a function of p, when (j> ~ v is not 
a fixed point of internal symmetries. Reheating is "instanta- 
neous" for p > 10/3 if Af is smaller than the value indicated 
on the plot. 



FIG. 15: Reheat temperature as a function of the energy 
scale during infiation, for p — 3 (black), 4 (blue), 5 (green) 
and 10 (red), when = is not a fixed point of internal 
symmetries. See main text for details. 
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FIG. 16: Relation between v and M from CMB normalization, for p — 3 (black), 4 (blue), 5 (green) and 10 (red), when (j> = v 
is not a fixed point of internal symmetries. Only the values for which Tr > GeV and v < mpi are shown. For p = 4, 5 and 10, 
the dots indicate the values of v and M beyond which reheating is instantaneous. The almost horizontal lines that are roughly 
superimposed at the top of the plot show the upper bound p6p on v below which preheating occurs in less than one oscillation 
of the condensate, for each value of p. 



We have studied the linear stage of preheating in the class of small field inflation models, where the curvature of the 
inflaton potential is negative during inflation. Although this is one of the most common classes of inflationary models, 
preheating after small field inflation remained much less studied than preheating after chaotic and hybrid inflation. 
On the analytical side, the problem is complicated by the fact that closed form solutions for the background evolution 
of the inflaton condensate are in general not available, while standard approximation like WKB for the evolution of 
perturbations are not accurate. Nevertheless, we saw that a detailed analytical description of preheating in this class 
of models is possible. The analytical methods that we developed in this paper may also be applied to the study of 
preheating in more complicated models. 

We showed that preheating after small fleld inflation is usually dominated by the tachyonic ampliflcation of inflaton 
fluctuations in the intervals of time when the inflaton condensate rolls in the region where the curvature of its potential 
is negative. A peculiar feature of this process is that the inflaton fluctuations experience a succession of exponential 
growths and decreases, so we called it "tachyonic oscillation". The exponential decreases of the fluctuations in these 
intervals of time arise because the coefficient of the decreasing mode is much larger than the coefficient of the growing 
mode, by an amount that depends on the scale of the fluctuation. Despite this temporary exponential decrease of the 
fluctuations, the full process is very efficient and we showed that preheating completes typically after less than five 
oscillations of the inflaton condensate. The range of scales of the fluctuations amplifled in the course of preheating is 
very wide, extending from the Hubble scale to the curvature of the potential at the minimum, with a peak given by 
Eq. (IST)) when preheating completes in more than one oscillation of the condensate. 

When the condition (l36l) is satisfied, the first tachyonic instability is so efficient that preheating completes in less 
than one oscillation of the condensate. In that case, the spectrum of the inflaton fluctuations at the end of preheating 
is peaked around the Hubble scale, or even slightly outside the Hubble radius because of the residual amount of 
inflation during preheating. Density fluctuations at the Hubble scale may then lead to an abundant production of 
primordial black holes, see also JlQ], which can put constraints on small fleld inflation models, see e.g. [40| for a recent 
update on observational constraints on primodrial black holes. The large field fluctuations amplified by preheating 
lead also to the production of gravitational waves (GW), with a peak frequency today 



where Hp and pp are the Hubble rate and the energy density during preheating, and fc,/a is the charactersitic 



VI. CONCLUSION 




(116) 
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momentum amplified by preheating ^ . GW from preheating after chaotic and hybrid inflation tend to have a frequency 
today above 10'^ Hz, which is too high to be observable by high-sensitivity interferometric experiments, either because 

pI^^ > 10" GeV or because K/a > Hp [13]. GW from the non-perturbative decay of condensates different from the 
inflaton, in particular super-symmetric flat directions, are more promising in this respect although the GW amplitude 
may be suppressed in that cases [l3]- On the other hand, GW from preheating after small field inflation may fall 
naturally in the frequency range accessible by ground-based and even space-based interferometers when fc*/a ^ Hp 
and the energy scale of inflation is small enough, which is common in this class of models. We also saw that preheating 
after small field infiation may be followed by a long matter-dominated stage before the universe thermalizes, which 
would further redshift and dilute these GW. A detailed study of the production of primordial black holes and GW 
from preheating after small field inflation is currently under way (soj . 

Finally, although an analytical understanding of the linear stage of preheating is often very useful to study its 
cosmological consequences, lattice simulations are eventually necessary to further study the non-linear dynamics. 
This is a difficult task in general for preheating after small field infiation, because of the wide range of scales that 
appear in the problem. In such cases, lattice simulations can be performed only for a limited range of parameters 
and the analytical results that we derived in this paper may be used to extrapolate the results to other regions of the 
parameter space. 
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Appendix A: Evolution of the Inflaton Perturbations in the Crater of the Volcano 



In this appendix, we solve the mode equation ()45|) in the interval [—tm, tm] and we match the solution to the WKB 
solutions (|49ll52p . As discussed in Section IIVI we do this by looking for an approximate Schrodinger potential for 
which we can find exact analytical solutions. 

In order to find guidance on the form of a suitable Schrodinger potential, let us first consider the mode equation 
(HSI) in the model plj) with p = 2 As discussed in Section |lTl this model in itself does not lead to small field 
inflation. However, in that case, the background evolution of the inflaton condensate can be solved analytically, and 
we can then flnd the analytical expression of the Schrodinger potential appearing in the mode equation. Thus we 
start by considering the rescaled inflaton potential 

{I (Al) 
instead of (|44|) . The conservation of energy analogous to (|43| then gives 



(^2 

when ^ (pf,. This can be integrated to give 4> = \/2/ coshi where, as in Section ITVl we normalize the time coordinate 
such that t = when (j> is maximal. The effective Schrodinger potential in the mode equation then reads 

-V"{cbit))^l~^j- (A3) 
cosh (t) 



which corresponds to (|53p for fcmax = 1, Uo = & and /i = 1. Note that the relation (|54|) is satisfled for these values of 
the parameters. 

Exact analytical solutions of the Schrodinger equation with this potential are known, see e.g. (36} . This potential 
provides also a very good approximation to —V" in the interval [—tm,tm] for the model ([57)) . see Fig. [Bl In this 



Eq. I|116| l is valid when the universe becomes quickly radiation-dominated after preheating. Note also that, in models where gauge fields 
play an important role in the dynamics of preheating, extra peaks may appear at well-distinct frequencies in the final GW spectrum 15j. 
The solution for the modes in that case has been studied numerically in the second paper of [2711 in the context of preheating after 
hybrid infiation. 
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approximation, the mode equation (j45j) reads 



It is convenient to define the parameter 



Uo 



cosh (/it) 



v = 



2 

max 



fc2 



(A4) 



(A5) 



and to make the change of variable 



-/i t 



(A6) 



Note that ly is real for k < fcmax and pure imaginary otherwise. In terms of the function x ^1"^ (1 — x) ^^"^ Vk, Eq. (|A4 
reduces to the canonical equation for the hypergeometric functions F{a, b, c, x) (4l| with 



(A7) 



b = 




The general solution may then be written as 



Vk 



= Ak x"/^ (1 - x)'"'^ F (a, b,c,l-x) + Bk x'''^ (1 - a;)""/^ F {c- b,c- a,l ~ v,l - x) 



(A8) 
(A9) 

(AlO) 



where Ak and Bk are constant coefficients, to be determined by the matching conditions. 

We can now derive the extra condition (|54p between the parameters by requiring that (jA10[) reproduces the known 
behavior Vk cj) oi the zero-mode fc = 0. At t = 0, the inflaton condensate reaches its maximal value and (j) changes 
sign, so that (jAlOp for fc = should be an odd function oft, i.e. it should be anti-symmetric around x = 1/2. One can 
check that this is not possible for the second term in (|A10[) . so that Bk = for the zero- mode. In order for the first 
term to be anti-symmetric around x = 1/2, i<"(o,6, c, 1 — a;) must reduce to a first-order polynomial, which requires 
a = — 1 or 6 = — 1 for fc = 0. Only the second condition is possible and, using (jA5| lASP for fc = 0, this gives 



/i2 



(All) 



which is equivalent to Eq. (I54p for fi. The numerical value is fi ~ 2.21. 

Let us now consider the behaviour of (|A10p in the vicinity of t = —t„i and t — tm, where the matching has to be 
done. Around t = — .t — >■ 1 where F{a, b, c,l — x) — > 1. This gives 



Vk 



Ak e"^* + Bke-''' 



fit 



for t 



(A12) 



Around t — tm, x — > where F{a, 6, c, 1 — x) can be expanded into a power series with coefficients that involves the 
Euler Gamma function 4l|. This gives 



Vk 



(AkCk + BkDk) e 



where the coefficients 



Ck 



{Ak Dk + Bk Ck) 



for t — > tr 



Dk, = 



s\n{Tiv) r(l - v)T{a)T{b) 
cos [(a - 6) 7r/2] ^ ^ _ 



sin {ttv) 



(A13) 

(A14) 
(A15) 



have been simplified by using the reflection formula for the Gamma functions. 
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1. Transfer Matrix for k <^ fcmax 

Consider first the case k < /cmax- In that case, we have to match (|A12p with (|5ip aX t — —tm and (|A13p with ([52 
at i = f„j. Ehminating the coefficients A^. and B}~ gives the transfer matrix 



ak\ /-Cfee-^^-^-A"*-. -Due'''' 



a*: 



(A16) 



where Xk is defined in Eq. ([57]) . Note that the determinant of the transfer matrix (|A16P is equal to one. 

As discussed in Section ITVl the WKB approximation is strictly valid only for the modes with k <C /cmax, which 
dominate the tachyonic effect. In this limit, the transfer matrix simplifies significantly. First of all, in the expression 
for CLk, the term in Ck 6^^*= at is neglibile with respect to the term in Dk e^'' bk for these modes. Next, the condition 
(jAlip means that b ^ ~l for k 0, so that T{b) oo and Cfe — >■ for these modes, see Eq. (jAl4p . From (jAl6p . 
this implies that \ak/bk\ ~ e^'^'' for A: — > 0, so that the solution ([5^ decreases exponentially during the whole interval 
of time [t,n,to]. Indeed, this was precisely how the condition (jAllI) was obtained, by imposing that the solution Vk{t) 
for k — 0, which increases exponentially in the whole interval [—to, —tm], be antisymmetric around t = 0. Now we 
consider the case < fc <C fcmax and compute the first order correction in k/kmax- In this limit 6 ~ — 1 — fc^/ (2/zfcmax) 
and we can use 

r(-l - e) = ^i— 4 ^ - for e^O (A17) 



to obtain 



e L^/c — -z — ; ; — -, ^ 



^ 1 (A18) 



2m fcmax sin {-Kv) T{1 - v) r(a) 

where the first factor can be calculated for fc = 0, which gives E ~ 2.15. To first order in k/k^^x, the factor Dk in 
(|A15p can also be calculated for fc = 0, which gives Dk — 1 after using (jAlip . Plugging all this into (|A16p leads to 
the transfer matrix (ISl 



2. Transfer Matrix for k ^ k„ 

For k > femax, we can use the formula above with 



u = iX , X = ^ (A19) 

M 

As discussed in Section ITVl the WKB approximation is valid for k ^ fcmax- In that case, we have to match (|A12p 
with (|49p at t — —tm and (|A13p with dSU]) aX t — tm- Eliminating the coefficients Ak and Bk gives the transfer 
matrix between the solutions and ([50)) for each oscillation of the condensate. The coefficients a'j, and /3^. in ([^^ 
after the j**^ complete oscillation of the condensate are further related to the coefficients cP^T^ and in ([SO)) during 
the oscillation by a-^ = a.-'fT^ e~*®* and = /3^^ e'^'', where 9^ is given by ([60]) with replaced by tm in the 
lower limit of the integral. We can then relate the Bogolyubov coefficients after and before the j"^ oscillation of the 
inflaton condensate as 




(A20) 



where 



^k = A/itm + / iOk{t)dt (A21) 



For k > fcmax, we take tm instead of ti in the lower limit of the time integrals in I I50II . 
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is the total effective phase accumulated during half a period of the condensate oscillation. The determinant of the 
transfer matrix (IA20p is again equal to one. 

The occupation numbers after and before the j"^ oscillation of condensate, nj, = and njT^ — can then 

be related by using (jA20p and the normalization condition — l/^fc~^P — 1- This gives 

nl = \D,\^ + {\C,\' + \D,\') rvir' + 2 \a\ \D,\ V^^f' ("f + 1) cos^i (A22) 

where 9^ + argCfc - ar gUfc + arg ^^.~^ ~ a rgaj.~ \ To calculate |£>fcp, we can use (|Alip and (|A19|) in (|A15p 

to obtain (p7|) . Finally, using (|A11I) and (jA19p in (jA14l) . we can simplify |Cfcp by using the reflection formula for the 
Gamma function and the fact that the complex conjugate of T{z) is given by T{z*). This gives |Cfcp = 1 + |-Dfcp and 
Eq. (|X22|) then reduces to Eq. (|66| . 
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